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Abstract 

Discrete and continuous non-intersecting random processes have 
given rise to critical "infinite dimensional diffusions" , like the Airy 
process, the Pearcey process and variations thereof. It has been known 
that domino tilings of very large Aztec diamonds lead macroscopically 
to a disordered region within an inscribed ellipse (arctic circle in the 
homogeneous case), and a regular brick-like region outside the ellipse. 
The fluctuations near the ellipse, appropriately magnified and away 
from the boundary of the Aztec diamond, form an Airy process, run 
with time tangential to the boundary. 

This paper investigates the domino tiling of two overlapping Aztec 
diamonds; this situation also leads to non-intersecting random walks 
and an induced point process; this process is shown to be determinan- 
tal. In the large size limit, when the ovcirlap is such that the two arctic 
ellipses for the single Aztec diamonds merely touch, a new critical pro- 
cess will appear near the point of osculation (tacnode), which is run 
with a time in the direction of the common tangent to the ellipses: this 
is the tacnode process. It is also shown here that this tacnode process is 
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universal: it coincides with the one found in the context of two groups 
of non-intersecting random walks or also Brownian motions, meeting 
momentarily. 
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Discrete and continuous non-intersecting random processes have given 
rise to critical ^Hnfinite- dimensional diffusions", like the Airy process and 
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the Pearcey process and variations thereof. These problems have been widely 
studied during these last ten years. 

One of these problems has been the study of the behavior of n non- 
intersecting Brownian motions on M, when n — )• oo. When the starting and 
end points of these particles are pinned down (say, for t = and t = 1), the 
cloud of particles, when n — )• oo, sweeps out in space-time a certain region. 
The Brownian motions have been shown to fluctuate like the Airy process 
near the generic points of the boundary of that region; see [26 1 [20 l [22l [29l |8] . 
The Airy process is governed by an Airy kernel, a double integral of a ratio 
of the exponential of cubic polynomials. Similar results have been obtained 
for the boundary of the frozen region of domino tiling problems |27] and 
random 3D partitions |25J. 

The boundary of that space-time region or boundary of the frozen region 
may have some singularities. A singularity can be a cusp, which corresponds 
to a bifurcation of particles or two sets of particles merging in the Brownian 
case; or a cusp in the boundary of the frozen region due to some non- 
convexity of the boundary of the model. The fluctuations of the particles 
near the cusp or the statistical behavior of the tilings are described by the 
Pearcey kernel, which is a double integral, with integrand given by the ratio 
of the exponential of quartic polynomials; see |28l ESI El H] , among others. 

Another situation is the one where two sets of particles meet momentarily 
and then separate again. Locally, the boundary of the two sets of particles 
has a singularity, which looks like two circles touching; such a singularity 
is a tacnode. It was an intriguing open problem to understand the local 
fluctuations of the particles in the neighborhood of this tacnode. Adler- 
Ferrari-van Moerbeke made an attempt in [3] by considering two sets of 
Brownian particles leaving from and forced to two points, and by letting 
first the number of one set go to oo, while leaving the second finite. This led 
to a rational perturbation of the Airy kernel at the point of encounter of the 
two sets of particles. Then, letting the number of particles in the other set 
of particles go to oo as well, one would expect to find the tacnode statistics. 
This approach remains an open problem! 

In [1], Adler- Ferrari- van Moerbeke resolved the tacnode problem for 
two groups of non-intersecting random walks {discrete space and continuous 
time); an explicit kernel was found, which is representable as a double inte- 
gral, for which the limit could be taken, thus leading to a kernel expressible 
-roughly speaking- as the sum of four Airy-like double integrals. 

Delvaux-Kuijlaars-Zhang |13j then found a kernel defined in terms of 
the solution of a 4 x 4 Riemann-Hilbert problem for two groups of non- 
intersecting Brownian motions (continuous space and continuous time). K. 
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Johansson [23] then found an exphcit kernel for non-intersecting Brownian 
motions, which seemed quite different from the one in 

The present paper deals with the random domino tilings of two overlap- 
ping Aztec diamonds: it investigates the fluctuations of the domino tilings 
near the region of overlap. The problem of a single Aztec diamond has been 
widely investigated by the combinatorics and probability community; the 
highlight was the existence of an inscribed arctic circle: inside the circle the 
domino's display a disordered pattern and outside a regular brick wall pat- 
tern; see [I51II61II81II71E21II9]. Johansson l22] then showed that the domino 
tilings near the arctic circle fluctuate like the Airy process, upon observing 
the boundary with an appropriate magnifying glass; the process is run with 
a time which is tangential to the boundary. This was done by showing that 
domino tilings of Aztec diamonds can be translated into non-intersecting 
random walks and a point process, which turn out to be a determinantal 
process. 

The problem of two overlapping Aztec diamonds is new. It translates 
into two distinct groups of non-intersecting random walks, which are also 
determinantal {discrete space and discrete time). In the limit, when the 
squares of the diamonds get smaller, we give three equivalent kernels for the 
fluctuations of the domino tilings near the overlap: the two first ones are 
closely related; the third one coincides with the one obtained for the two 
groups of non-intersecting random walks above, as in [4j; similar methods 
will be used to obtain the result. The latter is very different from the first two 
ones. Surprisingly so, they differ due to the use of the Christoffel-Darboux 
formula in the proof for one representation and not for the other. The first 
one is shown to coincide with the one found in |23| for the two groups of 
non-colliding Brownian motions. 

This paper is a first step towards a universality result for the statistical 
fluctuations near a tacnode, let it be a tacnode in the boundary of a frozen 
region or a tacnode for two groups of random processes meeting momentar- 
ily, including all variations on that theme, discrete and continuous time and 
discrete and continuous space. Finally, A. Borodin and M. Duits |9] investi- 
gate a Markov process with interlacing particles, which leads to a tacnode, 
but which is a different situation from our paper. 
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1 Double Aztec Diamonds and main results 




Diamond B 



Figure 1: Double Aztec diamond of type [n, m) = {7, 2) witli #{iiiIierH} = 
M = 2771 + 1 = 5. 

Domino tilings of double Aztec diamonds. Consider two Aztec dia- 
monds (checker boards) A and B of equal sizes n (= number of squares on 
the upper-left side) with opposite orientations; i.e., the upper-left square for 
diamond A is black and is white for diamond B. The diamonds A and B 
overlap a certain amount, according to the pattern indicated in Figure 1. 
Let M =: 2m + 1 (taken to be odd) be the number of white squares, belong- 
ing to the left-lower boundary of the double diamond A\JB; they correspond 
to the labels ii, . . . ,i2m+i in Figure 1. M also equals the number of black 
squares on the upper-right edge ol A[JB. The amount of overlap is given by 
n — (2m + 1). Cover this double diamond randomly by domino's, horizontal 
and vertical ones, as in Figure 2. The position of a domino on the Aztec 
diamond corresponds to four different patterns, given in Figure 4 below: 
North, South, East and West. 

Together with this arbitrary domino-tiling of the double Aztec diamond 
A\J B, one defines a height function h specified by the heights, prescribed 
on the single domino's according to figure 4 below. Let the upper-most edge 
of the double diamond Au B have height /i = 0. Then, regardless of the 
covering by domino's, the height function along the boundary of the double 
diamond will always be as indicated in Figure 2, with height h = 2n along 
the lower-most edge of the double diamond. Away from the boundary the 
height function will, of course, depend on the tiling; the associated heights 
are given in Figure 2. 
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The height function h obtained in this way defines the domino tihng in 
a unique way, because a white square together with its height specifies in a 
unique way to which domino it belongs to: North, South, East and West; the 
same holds for black squares. For example, a white square with a constant 
height along the four edges can only be covered by a North domino. 

The dual height function is obtained by polar refiection through the 
center of the four domino's, leading to the interchanges NorthoSouth and 
East -f-)- West; compare Figures 4 and 5. The dual height function leads to 
different boundary conditions on the double Aztec diamond, as shown in 
Figure 3. 



height and level lines for n— 7 



I 
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Fig 2. Domino-tiling of a double Aztec-diamond, the associated height 
function h and its level lines. 

Fig 3. Domino-tiling of a double Aztec-diamond, the dual height h and 
its level lines. 

Extension of the tiling region. It is convenient to extend the upper-left 
side of the Aztec diamond A by means of horizontal South dominoes and 
doing the same to the right of the diamond B, as depicted in figure 8. To 
be precise, at the upper level, one adds n South domino's, at the next level 
71 — 1, etc... So, the extension has a fixed tiling. Then, extending the height 
function to this new region, one has that the whole upper edge of the double 
Aztec diamond has height = and the whole lower edge height = 2n. 
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Figure 4. Height function h on domino's and level line. 




West 



Figure 5. Dual height fimction li on domino's and level line. 



The 2n outlier paths: denote by Ok the 2n level lines for this height func- 
tion at levels h = k — 1/2 for 1 < A; < 2n. These paths will be called the 
^^outliers'\ Since these paths correspond to distinct levels, they do not in- 
tersect each other. They consist of two groups, an upper-group of n paths 
oi,...,On, starting from the middle points of the left-edges of the black 
squares on the the upper-left side of domino A and ending up at the corre- 
sponding points to the right of the same diamond A, as in Figure 2. Then 
there is a lower-group of paths o„+i, . . . , 02n starting from the middle points 
of the left-edges of the black squares on the the lower-left side of domino B 
and ending up at the corresponding points to the right of the same diamond 
B, also as indicated in Figure 2. 

Given the extension of the tiling region, the outlier paths oi, . . . , are 
now extended by horizontal lines on the left and the paths On+i, . . . , 02n by 
horizontal lines to the right; they remain level lines for the height function. 

These outlier paths will now define two kind of particles, blue "dot- 
particles" and blue "circle-particles", according to the recipe indicated on 
the left hand side of Figure 6: put a dot-particle on the left-hand side of the 
segment of path traversing West and South domino's and a circle-particle 
on the right hand side of that segment, possibly with superposition of dots- 
and circle-particles, when edges of dominoes are in common; see Figure 8. 
Notice that a North domino will never be visited. Both the upper-outliers 
and lower-outliers thus have 

^{circles} = ^^{dots} = n + 1. 

The 2m + 1 inlier paths i-m, ■ ■ ■ ,im are the non-intersecting level lines 
for the dual height function h corresponding to the levels h = k — 1/2, for 
1 < A; < 2m + 1. This height function h is prescribed for single domino's in 
Figure 5, and is applied to an arbitrary domino tiling of the double diamond 
AL) B. Setting the height /i = at the lower edge of AL) B, the heights of 
the boundary edge is also completely determined, whatever be the tiling; 



7 



South 



Figure 6: Prescription of the "outher" paths, the (bhie) dot- and 
circle- particles on single domino's, and the map to the corresponding random 
walk 








C3 





Figure 7: Prescription of the "inher" paths, the (red) dot- and circle- 
particles on single domino's, and the map to the corresponding random walk 



the height of the upper-edge of ^ U i? is then h = 2m + 1, as indicated in 
Figure 3. The inher paths depart from the middle points ii, . . . ,i2m+i of 
the edges of the 2m + 1 white squares on the left boundary of the diamonds 
A and end up at the corresponding locations on the diamond B. 

Here also, one defines two kind of particles, red dot-particles'^ and red 
^^circle-particles" , according to the recipe indicated on the left of Figure 7: 
put a circle-particle on the left-hand side of the segment of path traversing 
East and North domino's and a dot-particle on the right hand side of that 
segment. The paths belonging to a vertical domino of type West contain 
neither circles nor dots. Here South domino's will never be visited. 

Circle- and dot-particles can then be superimposed when edges of domi- 
noes are in common. As boundary condition, one always puts a dot at the 
initial points of the paths «i, . . . , i2m+i and a circle at the end points. So, 
as depicted in Figure 9, each path has therefore 

^{circles} = ^^jdots} = n + 1. 

The inlier paths will not be used in the statement of the main Theorems, 
but they will play a crucial role in the arguments, which will first be made 
for the inlier paths and then, by duality, for the outlier paths. 
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Remark The height functions h and h are not the standard ones, as 
exhibited in [15]. They can be obtained from one another by an affine 
transformation. 

The system of coordinates {£,Y), the dot and circle particles are 

defined as follows: consider a horizontal "time-axis" i and oblique axes Y^'s 
through the values < £ < 2n + 1; namely, (i) axes Y2r-i = {Br,Dr), for 
1 < r < n, traversing white squares and (ii) axes = (vlr+i, Cr+i) for 
1 < r < n traversing black squares. The axes Y2r (resp. l2r-i) take on the 
values —n — m < k < n + m given by the location of the dots (resp. circles), 
obtained by moving the dots (resp. circles) belonging to the line (halfway 
in between the lines 1^2r-i and Y2r) horizontally towards the axis Y2r (resp. 
towards the axis Y2r-i)- The value of the dots and circles is specified by the 
equidistant horizontal lines beginning with the first one passing through the 
middle of the bottom most domino (level Y = —n — m) and the highest one 
passing through the middle of the top most domino (level Y = m + n). 

In other terms, an integer corresponds to a gap on the l2r-i-line, when 
the l2r-i-axis traverses a domino, not containing a circle-particle along that 
axis; similarly, an integer corresponds to a gap on the l2r-axis, when the 
y2r--axis traverses a domino, not containing a dot-particle along that axis. 
See Figure 8 for an example, involving the domino tiling given in Figure 2. 
So, the locations of the gaps on the axis l2r-i) say, indicate a change in 
pattern precisely at those locations. 

Another way of describing the dot particles on l2r and the circles on 
Y2r-i is as follows: put a dot on Y2r, in the middle of the black square, 
whenever the axis Y2r intersects an outlier path in that black square. Put a 
circle on l2r-ii in the middle of the white square, whenever the axis l2r-i 
intersects an outlier path in that white square. 

A weight on domino's, a probability on domino tilings and on 
non-intersecting random walks: put the weight < a < 1 on vertical 
dominoes and the weight 1 on horizontal dominoes, so that the probability 
of a tiling configuration T can be expressed as 

^#vertical domino's in T 

P(domino tiling T) = ^. . . — — (1) 

^ ^#vertical domino s m r 

all possible tilings T 

The definition ([T]) of the probability measure on domino tilings is not 
affected by the extension of the tiling region, since the extension has a fixed 
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Figure 8. Outlier non-intersecting domino paths, circles and dots. The coordinate 
Y2r (resp. Y2r-i) records the location of the dots (resp. circles), obtained by moving the 
dots (resp. circles) in between Y2r and Y2r-i to the axis ¥2,- (resp. Y2r-i)- The values 
thus obtained correspond exactly to the location of the dots (resp. circles) on Y2r (resp. 
127 -1) in the lattice paths description of Figure 16. E.g. for r = 5, Y2r takes on the values 
—9 — — rt — m<k<n + m = 9, except for gaps at —4, — 2, 0, 4, 6. 




Figure 9. Inlier non-intersecting domino paths ii, . . . , igm+i- Here the coordinate Y2r 
(resp. Y2, — 1) records the location of the red dots (resp. red circles), obtained by moving 
horizontally the red dots (resp. red circles) in between Y2r and y2r-t-i to the axis Y2r 
(resp. y^r+i). The values thus obtained correspond exactly to the location of the dots 
(resp. circles) on l2r (resp. l2r+i) in the lattice paths description of Figure 15. E.g. Y2r 
takes on the values —4, — 2, 0, 4, 6. 
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tiling. This domino measure induces a probability on the outlier configura- 
tions. This point process is a determinantal process, whose kernel will be 
given in Theorem This measure also maps to a probability measure on 
non- intersecting random walk paths. The duality between the "inliers" and 
"outliers", to be shown in Section [2j will enable us to obtain the formula of 
the correlation kernel for the outlier point process. 

The Frozen Region for the Double Aztec Diamond, when n — > oo. 

It is known from Jockush, Propp and Shor [18] that in the limit for a single 
diamond the stochastic region lies inside an "arctic" ellipse and the frozen 
region outside that region; if, say, diamonds A or B would be by themselves, 
the ellipses would be given by 

\ = 1, with q = and p = 1 — q = t, 

p q a + a ^ a -\- a ^ 

in terms of the weight a and in terms of the coordinates centered in 

the middle of the double diamond AL) B rescaled by 1/n. Thus the frozen 
region for each individual diamond would lie outside these ellipses. If there 
would be not much overlapping of the double Aztec diamonds AuB, so that 
the two ellipses have no point in common, it seems reasonable to assume 
that in the limit the frozen region will be outside these two same ellipses. 
Augmenting the amount of overlap will bring these two individual ellipses to 
a point where they merely touch, as indicated in Figure 10. Geometrically, 
this would be the case when the amount of overlap n — 2m — 1 of the two 
diamonds, rescaled by 1/n for n large, equals 

n-2m-l 2m + 1 2 ^ 

= 1 ~ 1 = 1 — 2vpg =: 1 — r. 

n n a + a ^ 

Indeed, the centers {'^m, ±{2m + 1)) of the diamond A and B, rescaled by 
1/n, are then given by 




a + a^' a + a ^ 

and thus these two ellipses are tangent to each other at (0, 0). So, one expects 
the frozen region for a domino covering of the double Aztec diamond Au B 
to be given by the region outside the two ellipses within AU B; this is not 
shown in this paper. Still this argument suggest the scaling 

m 1 



n a + a ^ 
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One might conjecture that the two regions, in between the elhpses and the 
boundary of the double diamond (having the origin in common), are frozen 
with North domino's. 




Figure 10: The macroHcopie description of the double Aztec diamonds 
and their frozen region outside the elhpses. 



Main results: In view of Theorem |l.l[ define the kernel for the one- Aztec 
diamond obtained from |22j by setting u i— t- —u, v i— )■ —v, 



j^oncAztcc(^2r,x;2s,y) 



(27ri) 



du 



7r3 



v-u u^-y (1 + avY-''{l - 



(2) 

with the integral ip2{s-r): defined below. 7^ refers to circles of radius r 
around 0. In Q and throughout the paper many different radii will be 
considered, subjected to the following inequalities: 

< a < r3 < r2 < S2 < si < Ti < S3 < and < a < p < 1. (3) 

The contours Fq or Fo.a refers to a contour containing or 0, a, but no other 
pole of the integrand; setting n even, we now define the following functions. 
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used throughout the paper: 

I)''-'' f , / dv (1 + at;)^(l - 



5(2...; 2s,y):J^j Jul 



(27ri)2 jT^^ J^^^u-v u^+^ ipa{2n;u) 
ly-y f , f dv (Pa{2n;v) 



u-v ny+'^+i (1 + auY{l - 



' (27ri)^ /^^ J^^_^ V - uv'^+^ ipa{2n,v) z - V 

.f^N = -/ |^(-.)V(2n;n), (n) = - / 

.- / -dz; (2) ._ f -dv Va{2n-v-^) 

^^^-fr,^2ni{v-z)M2n,v)' ''^ ^""^ - %,2.i{v-^-) {-vf^^ 

^ (4) 

For future use, g^'^\n) is a vector whose £th component is g^^\n); similarly, 

/i *■*•'(.) is a vector whose ^th component is h^^\.). Also set (sometimes p will 
be inside and sometimes outside) 

H('\z-^) := det(l - = det(/ - KW^^-^)) 

= det(l-K(i)(z-i))p, (5) 
with K^^"^ = X[p,p+i,...] k/X[p,p+l,...]- 
The same notation will be used for kernels on the reals, like 
det(l-i^), = det(l-X[,^^)ifXf,_^)) and (l-K)-' = (l-X[.,oo)i^^[a,oo))~'. 

We now state one of the two main Theorems, where P refers to the proba- 
bility ([1]), defined above. The forms of the same kernel obtained in Theorem 



1.1 are quite different; it is the use of the Christoffel-Darboux formula at 



some point of the proof which is responsible for the seemingly entirely dif- 
ferent forms (Isj) and ( 10 ) of that same kernel. The main event for which the 
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probability will be computed in this paper concerns the dot-particles, thus 
belonging the even lines Y2r- The following events are all equivalent: 

{The line Y2r has a gap D [k,£]} 

:= {The line l2r has no dot-particles along the interval [k,i] C 1^2r} 

{The domino's covering black squares along the | 
interval [k,i] C l2r are oriented left or down J 
= {The height function h is flat along the interval [k,£] C l2r} 

{The line does not intersect any of the outlier 1 
paths in the black squares along [k,£] C Y2r J 



The line 1^21- has a gap D {k} means: at that location one of these two configurations holds 



k+i 



k + l 
k 



The line has a gap Z) {k. k + 1}: at that location one of these three conhgmations holds 

Figure 11: The line has gaps D {k} or {k, k + 1}. 



In Theorem |1.1[ one shows that the dot-particles along the lines Y2r 
with even indices form a determinantal process. The circle-particles and 
the combination of circle- and dot-particles form, at odd and even line Yi 
respectively, a determinantal process with a kernel which could be computed 
as well. 



Theorem 1.1 The dot-particles for the outlier paths form a determinantal 
point process with correlation kernel given by ^) under the probability 

measure induced by In particular, given integers —n — m < k < £ < 
n + m (n even), the following probability holds for the single line Yn (Y -axis 
in the middle) : 



' {The line Yn has a gap D [k,£]) = det I 1 — X[k,e\ 



T.X[k,t] 



(6) 



14 



where Kn^mix,y) := W^^{n, x]n,y), the extended kernel given below. For 
multiple lines ¥2^, 1 <i < s, 

F {the line has a gap D [/sj,^^]}^ 

= det I 1 - X[k,A]K^lni'^ri, xf, 2rj,Xj)x[k,/A 

K^^^ is a perturbation of the one-Aztec diamond extended kernel, with an 
inner-product involving the resolvent of the kerne^K, defined in 

(-l)^-'^K-^(2r,x;2s,y) 

= ^'""(2(n - r + l),m - X + l;2{n - s + l),m - y + I) (8) 

+ ((1 - K)2i+ia^y,s{k), b^xAf^))e2^2m+i,...) ' 

with 

j^oneA.tec(2(„ - T + 1 ) , m - rc + 1 ; 2 (u - s + 1 ) , 771 - y + 1 ) 

= -ls<r{-ir-'^i^2{s-r){x,y) + S{2r,x;2s,y) 

An alternative expression is given by: 

i-ir^y^^p±^K^^l^2r,x-,2s,y) = -l.o.(-ir-V2(.-.)(^, 

Sxj^y f dz 



+ 



+ 



27ri ./ro 



/ rf^ / ^ (1 _ rW{z-')){1 - 



(10) 

\ / ^2(2m+l,...) ^^^^ 



^{.f{^)ygif^))e'^{2m+i....) = X]2m+i fi^)9{^)- Also remember the notations in Jsk 
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Given the weight < a < 1 on vertical dominoes, the following quantities 
vq < 0, and A, p, 9 > will come into play: 

1-a ,o a(l + a)5 ^ {1 - a^ f/^ 
vo ■= --— < 0, := ' p := -Avo = f— > 

61 := V/9 a + a-i = > 0, and so = 61- — „. . 

V a / (1 — aj(l + a^j 

(12) 



Consider the scaling ( 13 ) of the following quantites: the type (n, m) of the 
double Aztec diamond, the space and time variables x and s in terms of the 
size n = 2t of each diamond, where a is a parameter, which measures the 
pressure between the two diamonds: (for future use set a := 

n = 2t, m= — + apt^/'^ 

a + a ^ 



X 

r 



2aHnt"^ + i^pt^'\ y = 2a'eT2t^/^ + ^^pt^'^ ^^^^ 
t + {l + a'^)eTit^'^, s = t + {l + a2)^r2^2/^ 



In the second main Theorem, three formulas will be given for the limiting 
kernel. To state the first two ones, the following functions are needed. Given 
the parameter s, define an extension Ai^*^(x) of the usual Airy function 
Ai(x); further an Airy- type kernel, defined by means of Ai*-'^^ and finally a 



new function ^l^: 



Ai(^)(x) := — [ dz e^'/3+^'^-^^ = e'^+i''Ai(x + s^), 

/"OO 

K'^'~^\x,y) := Ai^"\x + u)Ai'^-^\y+u)du, with Km = i^i?'^^ 
Jo 

= / dn e-(^-")"Ai(x + + u)Ki{y + /^^ + 

AI{k) := Ai(^)(e + 2^/=^^) - / Ai(^)(-^ + 2^/^l3)k\{K + I3)d(3 

Jo 



a=l3=0 
u), 

(14) 



and 



p(t;Ci,6) := - 



(ei-e2)^ 



a/4 



ITT 
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Also define the extended kernel for the Airy process (see [22l [23] ) 

POO 

j^Ai.yP.ocess(^^^ ^ ) ^ / g A(r2 ) + A)Ai(e2 + X)dX 

Jo 



(15) 

The following expressions will be used in the second formula for the limiting 
Q(k) := [(1 - xaKAiXa)-\aAi] (k) 

POO 

Q(C) := y_ Q{K)e^'''''^dK (^g) 

Ja Jc 
together with the following functioiij^ 







Ai(2-^''^«)(K + 2-1/3^) 



/■OO 

e(s,0:=2"^/M (iKQ(K) 

/■OO roo 

+ / dAQ(A) / daAi(a + A)Ai(2"'^'^)(a + K + 2-1/3^) 



(17) 



The tacnode process measures the probability that the height function 
h is flat along small intervals of length 0{n~'^/^) in the y-direction (see 
Figure 8), which themselves are 0{n~^^^)-close to the center of Figure 8; 
the tacnode process is run with a time scale of the order 0{n~^/^), near the 
halfway point i = n (see Figure 8); i.e., near the tangency point in Figure 
10. 



Theorem 1.2 Given the scaling (13), depending on the arbitrary parameter 
a, with a := > 0, the scaling limit of the dot-particles, along the lines 

Y2r, near the tangency point of the two ellipses is a determinantal process 
with kernel: 

hm i-voy-^^'-'i-iy^K'^^W^r, x; 2s, y)pt^/^ = K*-(ri, 6; T2, 6) (18) 

E— >-00 

where the tacnode kernel 'K}^'^{ti,^i;t2,^,2) has three equivalent forms. 



^(^ -f-)- — ^)means; replace ^ — >■ — ^ in the prior expression. 
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(i) A first form is a perturbation of the Airy process kernel (15): 



+ 2^/3 ((1 - Km)^Uii^^) iX)A-;i^{X)dX. 

(19) 

where E] 

-lr,>r,p{n - r2;6,6) + K'^r"'\a - 6,^ - 6) 

(20) 

(ii) It can also be expressed as 

]K*^^(n,ei;T2,6) 
= -lri>T2P(n - T2;6,6) 

+ 25^" [((l-KAi)^U^i_.) (A)Ai(--^)(e2 -a + 2^/3 A) + {e. o 

(21) 

(iii) Another equivalent limiting kernel is given by the sum of the four double 
integrals, where 6 > is arbitrary: 

I^*"'(ri,e2;r2,6) 

= - l[r2<ri]P(n - T2;6,6) + ^(ti - T2,Cl - 6) 



dX. 



(27ri)2 y^+ijj gi^-.. Ve«2" e-?^- 



3 

, , eV-'^^ e^i"Ve^'" e-«i"\ (1 - 5'(ii))Q(-?;) 



(27ri)2 y^+iR g-^-<x^; e-^'^' Ve^^'' e-«2^y u-v 



„3 2 



1 f ^ f ^ er^^e^/e^ e-«l"^^ (1 - a'(-v))Q(u) 



(2vri)2 y-25+iR pf-a^ e-2^He«''' e'^^V ^ " 



e 3 

"■^ ,,2 



1 /",/■, e---'^"e^i" /e«i" e'^i^X Q(u)Q(-w) 
(27ri)2 y^+ij, g-^^-a. \e^^- e-^^- ) u-v 



(22) 



Consider two groups of Brownian motions, leaving from and forced to 
distinct points at times t = and 1, tuned in sucli a way that they meet 



Vith q,y{r, £,) := e'^ 



(<^-«) + (2/3)t 
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momentarily at t = 1/2. Then, when the number of particles gets very large, 
so as to preserve this tuning, the statistical behavior of the particles near 
the point of encounter are governed by the following kernel, as is shown by 
Johansson in [23] 



+ 23 ((l-KAi)^^yi:g_,) (k) [A\_^iK) - Ai(-i)(-ei -a + 25^) dK 

(23) 



Theorem 1.3 The Brownian motion tacnode kernel (23) is equivalent to 
anyone of the kernels in Theorem 1.2, 

It follows from the above result that the correlation kernels obtained by 
Adler, Ferrari and van Moerbeke in [4], by Delvaux, Kuijlaars and Zhang 
in [13] and by Johansson in |23j are all equivalent, i.e., they all describe the 
same tacnode process! 



2 From Aztec Diamond paths to Lattice paths, via 
Zig-zag paths 

Duality and Zig-zag paths. The two sets of lattice paths are dual in 
the following precise sense. The circle-particles (dot-particles) for one set 
of lattice paths fill the circle-gaps (dot-gaps) for the other set; see Figure 
17. To see this, we define the zig-zag paths, initiated in [ISl UHl [20],. There 
are two types, zig-zag paths around black squares and zig-zag paths around 
white squares. 

Remember the axis Y2r, connecting C^+i and A^+i, as in Figure 8, tra- 
verses 2n + 2m + 1 black squares of the double Aztec diamond; it does so 
in two different ways, which can be described as follows: the black square 
can belong to a South or West domino, in which case the height h goes up 
by 1 from top to bottom, as in the left hand side of figure 12 below; put a 

*The kernel below differs from Johansson's kernel in (23] by a transpose and a conju- 
gation by an exponential; i.e., multiplication of the kernel by . 
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blue dot-particle to the left of the oblique line. Or the black square belongs 
to a North or East domino, in which case the height of the oblique line 
remains unchanged, as in the right hand side of figure 12 below; put a red 
dot-particle to the right of the oblique line. 



W or S N or E 




Figure 12. Zigzag paths around black squares 

So, in the left configuration, the line, running from top to bottom goes up 
by /i = 1 and for the second, it stays flat. Since the line must go up from 
h = to h = 2n, the line must traverse 2n black squares of the first type, to 
reach height h = 2n and must stay flat along the 2m + 1 remaining squares, 
since the total number of black squares the oblique line traverses equals 
2n + 2m + 1. Therefore the red dots appear there, where the blue dots are 
absent. Moreover the blue dots must belong to W or S domino's and the 
red dots to E and N domino's. 



W or S N or E 




Figure 13. Zigzag paths around white squares 

Similarly, the axis l2r-i connecting the points Dj. and (as in Figure 
8), traverses 2n -\- 2m + 1 white squares, also starting at height h = until 
it hits height h = 2n. In order to reach that height it must traverse 2n 
white squares of the left type of flgure 13, for which the oblique line goes 
up by /i = 1 and stay fiat for the others (white squares of the right type 
in Figure 13); so, these two situations are dual to each other. Then put a 
blue circle-particle to the right side of the left square of Figure 13 and a red 
circle-particle to the left side of the right square of Figure 13. The white 
squares, carrying the blue circle-particle (left of Figure 13) must belong to 
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W or S domino's, and those carrying a red circle (right of figure 13) must 
belong to N or E domino's. 

So the axis Y2r, as in Figure 8, records the location of the blue dot- 
particles of the line [C^+i, and the axis Y2r-i records the location of 
the blue circle-particles of the line [Dr , Br] ■ To show that the blue dot- 
and circle-particles belong to the level paths for the height function h, as 
constructed in section 1, it suffices to notice that the blue circle-particles 
belong to S and W domino's and the blue dot-particles as well, from the 
previous considerations. Upon connecting them, they are part of level paths 
for the height function h. This shows that the two constructions, the one 
using the zig-zag paths and the one using the arguments of section 1 are the 
same. One does the same for the red dot- and circle-particles, which have 
been shown to be dual to the blue particles. 

Consider now the red dots and red circles on the line halfway in between 
and Y2r+i, as in Figure 9. Move the red dots horizontally to the left 
towards the axis Y2r and the red circles to the right towards the axis l2r+i- 
Then in comparison with Figure 8, the axes l2r and i2r+i will have red dots 
and red circles precisely at the sites which have no blue dots and circles, 
showing the duality. 

Outlier (blue) circle-particles 




S W N E 



Outlier (blue) dot-particlcs 




S W N E 



Outlier combined (blue) eirele-particlcs and dot-particles 




S W N E 



Figure 14a. The three paths above are level paths for the height function 
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h, as in Figure 4. 

Outlier (red) circlc-particlcs 




S W N E 



Inlicr (red) dot-particlcs 




S W N E 



Outlier coiirbincd (red) circle-particles and dot-particles 




S W N E 



Figure 14b. The three paths above are level paths for the dual height 
function h, as in Figure 5. 

The extended outlier domino-paths are transformed to non- 
intersecting paths on a lattice, by means of the correspondence of Figure 
6. Here each path on the double Aztec diamond is made of segments covering 
West, South and East domino's. The path on the lattice is obtained by 
sewing together the corresponding segments on the right hand side of Figure 
6, with the accompanying dots and circles. To summarize, the outlier lattice- 
paths will thus contain 2n + 1 steps starting with a dots and ending with a 
with a circle, so that each path has #{dots} = ^{circles} = n + 1. 

Notice that the location of the circles on the axis l2r-i and the dots on 
in the double Aztec diamond picture (Figure 8) corresponds exactly to 
the circles and dots on the axes l2r-i and l^r of Figure 16. The map of 
Figure 6 from domino paths to lattice paths was designed to achieve this 
precise correspondence between the Y^'s in the domino description and the 
lattice description. 

The inlier domino-paths are transformed to non-intersecting 
lattice paths, by means of the recipe in Figure 7. Indeed, each path 
h, ■ ■ ■ , ^2m+i on AL) B consists of segments covering East, North and West 
domino's, together with their appropriate dots and circles; the recipe is to 
map the paths on the domino to a path on a lattice by sewing together the 
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Figure 17: Lattice paths: superimposing in- and outliers. 

paths on the right hand side of Figure 7; this leads to a path on the lattice 
with alternating circles and dots. As is suggested by Figure 7, a path on a 
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West domino gets transformed in a vertical line (of unit length) on the lat- 
tice, possibly repeated as many times as the path traverses West dominoes. 
The dots as initial boundary condition for the inlier paths translate into a 
dot with an added horizontal segment on the lattice. The final boundary 
condition for the inlier path, namely the circle translates into simply putting 
a circle at the end of the path on the lattice. 

One thus distinguishes two kinds of steps for the inlier lattice paths: 

A-steps = {steps from o to •} = {steps from T-xorx— t-x + I} 
5-steps = {steps from • to o } = {steps from x ^ x — k, with A; > 0}. 

3 Expressing the kernel in terms of orthogonal 
polynomials on the circle 

Consider the following weights on the unit circle in C, 



and 



Pai^)=(Y^) ^^dp„«(z) = ^^— (24) 
Paiz) := p!^{z)p^{z) = (25) 

\^ z) 

The following operation * = * will be used throughoulj^ 

f{x) * g{x) := f{x) * g{x) := ^ f{x)g{x). (26) 

Also remember from Q the definition for < s G Z, and x,y G Z, 

which satisfies the semi-group property with respect to the operation *, 

'ilJ2s{x,-)*'^2r{-,y) = '4'2{s+r){x,y), for s, r G Z. (28) 

We now state the Proposition below, which uses arguments similar to the 
ones in 0]. 



^The following formula involving the *-operation will often be used, e.g. in (28 1 
j / dz f dz 
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Proposition 3.1 The point process on the even lines associated with the 
two groups of non-intersecting outlier lattice paths, corresponding to the dou- 
ble Aztec diamond of type {n,m), is a determinantal point process with the 
following correlation kernel: 



(29) 



K.ZU2r,x-2s,y) 

= -ls<rV'2(s-r)(a;,y) + ^n-2r(a;, •) * K^'^mC"'' • ; o) * V'2s-n (° ,2/), 

expressed in terms of the kernel at the half-way axis 2r = 2s = n : 

^^n^ni^^X ■,n,y) = Kn,m{x,y) := (1 - Kn,m){x,y), 

where 

(30) 

with bi-orthonormal polynomials Pj. and Pk on the unit circle, for the weight 

Proof: Step 1: The kernel at time 2r = 2s = n. Remember the weight 
< a < 1 for vertical dominoes and the weight 1 for horizontal dominoes. 
We first consider the inlier lattice-paths, as depicted in the right hand side 
of Figure 7 and Figure 15. At the level of the inlier lattice-paths, the weight 
will figure in only when x — )• x + fc for k ^ 0. Therefore, by ^], the following 
weighting holds for j = 0, . . . , n for i?-steps (from dots to circles) and for 
j = 0, . . . , n — 1 for ^-steps (from circles to dots) 



5-steps : tp2j,2j+iix,y) ■-- 



a^~y, if y - X < 1 f dz z^^'^ 



otherwise jr-n 2tx\z 1 



A-steps : V2i+i,2i+2(a;,y) := < 



a, if y — X = 1 
1, if y — X = 
0, otherwise 



1, ify-x = 0)> = / + 

(31) 



rn 



Then define 

■fAns = V^r.r+l * ... * -05-1,5, if S > r 

= 0,ifs<r. 



'B-steps for outlier paths correspond to A-steps for inlier paths. 
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Therj^ for instance, for k even, one needs | -B-steps and | A-steps to go from 
—m + i — 1 to X in time k, where — n — m<x<n + m and 1 < i < 2m + 1, 

M-m + ^ - 1, x) = / A^,-m+^-l-^ (l±^^ (33) 



ro,ii 



2T:iz 



and n + 1 — I i3-steps and n + 1 — | ^d-steps to go from x to — m + i — 1 in 
time N — k: 



— + - 1) = ^ /^^^^ (1 - |)(2n-..2)/2 (34) 

Then from the equivalence between the double Aztec domino configurations 
and the inlier paths, from the fact that the weighting (31) follows from the 
probability ([T]) of a domino tiling and from the Lindstrom-Gessel-Viennot 
formula for non-intersecting paths, it follows that the probability (inherited 
from the step weights (31 )), of the paths being at xi, . . . , X2m+i at time r is 
given by 

F{r;xi, . . .,X2m+i) 

= ^r-^ det{'4)o,r{-m + i - ^, Xj))l<ij<2m+1 det{'llJr,2n+l{Xj, -m + i - l))l<jj<2m+l 

Hence, by the Eynard-Mehta Theorem \14:\ 120], the extended kernel is given 
by (Note lr-<<i can be omitted in view of definition (32) of ilJr,s) 

2m+l 

= ^ iJr,2n+iix,-m + i-l){[A'^]ij)-ipo^si-rn + j-l,y)-lr<:s-tpr,six,y), 

i,j=l 

(35) 

with the entries of the (2m + 1) x {2m + 1) matrix A, 

Aij = tpo,2n+i{-m + -m + j - I) = fi* gj, 

defined in terms of the functions 

/,(y) = V'0,n(-m + j-l,y)= <[ -^z-™+^--i-V^(z) 

r T (36) 

g,{x) = V„,2n+i(x, -m + i-l)= i -^z^^+i+Vf (z). 



'^The subscript will sometimes be omitted in V'o,fe, which makes it compatible with 
formula ([27 1. 
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with Pa{z) and Pa{z) as in (24). Consider now the kernel at a single (half- 
way) time n, which we already assumed to be even: 



i{x,y) := K^^jnii^,x;s,y) 



2m+l 

g,{x){[A-%)f,{y). (37) 



Then one performs row operations in the determinants of the Lindstrom- 
Gessel-Viennot formula; this amounts to taking linear combination fk{x) 
and gk{x) of the fi{x) and gi{x) with 1 < i < k, 



fkiy) 



gk[x) 



1 
27ri 

1 
27ri 



(38) 



dz 



p'^{z)z^^^Pk-i{z-') 



ro.a 



such that the new matrix A^i in (37) becomes the identity; i.e., 



fk{x) * gi{x) = ^ fk{x)gt{x) 



1 



dz 



Pk-l{z)Pe-l{z-^)pa{zy- 

ro,a ^ (39) 



thus leading to orthogonal polynomials on the circle for the inner product 
^ , » with regard to the weight pa{z)^, as defined in (25). With these 
new functions fk{y) and gk{x), the kernel (37) has the simple form: 

M 



Kn,m{x,y) = Ygk{x)fk{y) 



k=l 



1 



(27ri)2 



dz 



-y-m 



2m 



' i 0,a i 0,a,2 A;=0 

(40) 

where it is legitimate to include z in the contour of integration, since the 
ti;-part of the integrand is holomorphic m. w ^ 0, a. According to [12], the 
statistics of the non-intersecting outlier lattice paths is then described by its 
dual kernel, 

Kn,m(a;, y) = 6x,y - ^n,m{x] y). (41) 



This establishes formula (29) for 2r = 2s = n, with kernel ^n,m{x,y) as in 
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Step 2 : The extended kernel. Comparing the form of the extended 



kernel K^^J^(2r, x; 2s, y) in (35) and the form of the kernel 'Kn,m{x,y) = 
K^^^(n, x; n, y) in (37), and using the semi-group relations (28) on the ■0- 
functions, namely 



•^n-lrix, ■) * ll)n,2n+l{-, -m + i - 1) = 'tp2r,2n+l{x , -171 + i - 1) 
i^O,n{-'^+j - Ir) * i^2s-n{-,y) = V'0,2s(-"l + j - 

one checks that (these arguments appear in Lemma 5.2 of [4J) 
K^^(2r,x;2s,y) 

= -Ir<s1p2{s-r){x,y) +1pn-2r{x,-) * '^"n^ni^, ■ ; n, o) * V'2s-n(o, y). 

The dual extended kernel is then given by 
K^^(2r, x; 2s, y) = 6rs6,y - K^^(2r, x; 2s, y) 

= 6rsSxy + h-<s'>P{2s-2r){x, V) " i^n-2r{x, ■) * K^^mC'^J S °) * V'2s-n(-, 2/) 
= {Srs^xy + '^r<s1p(2s-2r){x,y) - 1pn-2r{x, •)!(•) °) * 02s-n(°,y)) 

+ tpn-2r{x, •)*(!- lKn,m)(-, °) * 1p2s-n{°, V) 
= -ls<r'<lj{2s^2r){x,y) + V'n-2r(x, •) * K„,m(- ; o) * 0(2s-n) (°, y) , 



(42) 



where one uses (41) and the identity 



5rs5xy + '^r<si'2s-2r{x, y) - ijn-2r{x, •) * o) * 1p2s-n{°, y) 

= -ls<T^{2s-2r){x,y), 



concluding the proof of Proposition 3.1 



Orthogonal polynomials on the circle and Fred- 
holm determinants 



In order to be able to use formula (|30|) for the kernel IK^ 

{x,y), one needs 



to express both, orthogonal polynomials and Darboux sums of orthogonal 
polynomials on the circle, in terms of Fredholm determinants of certain ker- 
nels. In the last part, we shall state some Fredholm determinants identities 
from 



28 



In order to do this, one first introduces the inner-product 



1 



t,s 



du 



f{u)g{u )e^J 



(43) 



27ri Jgi u 

upon setting t := (ti, • • •) ^ and s := (si, S2, . . .) £ The inner- 
product ^ , S>, as in (39), is a special instance of ( , )t,s upon picking the 
special locus 



£ = {(t, s) with tj 



-nt^ and 



(n+l)y}. 



Indeed, 



and thus 



(1 -han)" 



n+l 



Also notice that (see definition (H 



{l + auni-^r+' = M'^n;u). (46) 



Lemma 4.1 Given the biorthonormal polynomials Pk and Pi, 



1 

2TTi 



^ dz 
Pk{z)Pe{z'^)pa{z)— = 6k,e, 

ro,a ^ 



the following identities hold: 
Pk{z) = 

^P,{z)P,{w-') = 

j=0 



^fc(l_|)"+iffW(^-i) 
^iffc(0)i?fc+i(0) 

tt;-''(l + aii;)-"ijf^(tf) 
V^fe(0)i7fc+i(0j 



(l + a?/;)"(l-f) fffc+i(0) 



where 



H^f}\z-^) := det(/ - /r(^)(z-i)),2(fc,,+i,...) 
i7f)(z) :=det(/-E:(2)(^))^^^^^^^^^^^^^ 
//f (z-i,«;) :=det(/-if(3)(^-i^^))^^^^_^^^^^^^^^^ 



(44) 



(45) 



(47) 



(48) 
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with 

^i^fi^^^) for\z\ > \u\, \v\ 



(-l)k+e r r , 1 z-uipa(2n,u) , , 
^ dni dv^ ^ , (49) 



(27ri)2 7|„|=p<i 7|^,|=p-i>i v''+^v-uz-vipa{2n,v)' 

Kj^j{w) for Iw^^l > \u\, \v\ 

t_l\k+e J- J- , 1 /w-^ -u\ipa(2n,v-^) 

au (b av- 



(27ri)2 J|„|=p<i J|^,|=p-i>i v''+^v-u\w ^ -vJ ipai2n,u i) 

(50) 



K\,1{z ^,w) /or a < r2 = < S2 < Si < ri = |z| < a ^ 
._ (-1)''^^ / , / . u^^^ 1 ('u-^)(v-w^)(/'a(2n,u) 



(27ri)2 f^^^^'^ J^^^^^ v'^+'^ V - u {v - z){u - w) (pa{2n,v) ^^^"^ 

= <i,.+i(0) + (|; - i)4'i,(^'^)/^a(-), 

being rank 1 perturbations of K^^\o) , 
as kernels in k and I. In these formulas, 



gf\n) = -j^ ^(_^)Va(2n;tx) = -jf ^(_n)-^- Va(2n; n-^) 
(2)/ s / du (-uY f du i-uY^"'^ 

g} '{n) = ^ - rts 



(52) 



an 



Tro^_^2Tri(fa{2n,u ^) Jr^^^ 2TTi (pa{2n,u) 
h^\z-') :=- 



dv 1 
r(,_^ 27ri(f - z) (-f )*=+Va(2n, -u) 

oo ^ 



a=0 



^To^a in the first integral can (and will be) be realized as 7^ with a < |?;| = s < \z\. 
Also To in the second integral can be realized as 7^, with |u| = s < jw"^!- Finally Fcu, in 
the third integral can be realized as 7s, with |m| < s — \u\. 
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hf\w) :-- 



dv 



2TTi{v — w^^) (— f )^+^ 
w du 



ro,i 



2TTi(u — w) 



-u) ipa{2n,u) 



(53) 



OL=0 



The unperturbed kernel can then he expressed in terms of the g^^^ 's: 



.(1) 



(2) 



'(2), 



(54) 



a=0 



and therefore K^-^^^ {0) = K^'^\0). 



Proof: It was shown in [6j (see also the lecture notes |30) ) that the function^ 



pl\t,s;w 



-1^ 



w 



k Tfc(t- [Z-1],S) 

Y^Tfc(t,s)rfc+i(t,s) 

-k Tkit,S+[w]) 



(55) 



■\/Tfc(t,s)rfc+i(t,s) 



are bi-orthonormal polynomials with regard to the inner-product (43), i.e., 
{P^k\pf ^)t,s = 5k,i and that 

|:^-'(t.s;.)^-'(M;.-).(^)- -^'^;i-;;i;yH) 

In the formulae above, the r„(t, s) are 2-Toda r- functions and are defined 
as a Toeplitz determinant, which is also expressible as a Fredholm determi- 
nant of the kernel (58) below, using the Borodin-Okounkov identity [11]. We 
obtain two formulas, one obtained from the other by substitution n i— )• 1/u: 



rp(t, s) := det 
= det 



1 

2TTi 

1 

2TTi 



du 



SI U 

du 

SI u 



i<k,e 



l<k/<p 



Z{t, s) det (1 - K(t, s)),.„p^p+,^...„ , Z(t, s) := e" ^ '^'^ , 

(57) 



For a £ C, one defines \a] = 
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and the kernel K(t, s) is given by two seemingly different expressions, upon 
using the two different expressions for Tp above, related by t -f-)- —s: 



K(t,s)fc, 



1 



(2^i)2 

1 

(2^ 



du 

|m|=p<i Jr 



dv 



u 



1 



|j)|=p-i>i 



du 



\u\=p<l 



dv- 



u 



1 eE°life«^+s.«-^) 



|jj|=p-i>i 



yk+l y — n 



(58) 

This is valid insofar Yl'jLiitjU^ + SjU ^) is analytic in the annulus (p, p ^). 
In view of the expressions (55) and (|56|), involving shifts, we now compute. 



using ~ — ln(l — v/z) (for \v/z\ < 1) in the second formula 



1 1 _ « eE=li{ti«^+«i«-^) 



K(t — [z ^],s)k,i^ with |z| > \v\ 



du 



K(t- [z-i],s+H) 



dv- 



u 



k£ with < \u\, \v\ and |ii|, |f| < \z\ 



1 - ^J:T=i(ij^'+^i^-') 



{2mf 



du 



\u\=p<l 



dv- 



u 



|i,|=p-i>i 



,A:+1 



v-u{l- ^)(1 - ^) eEr=i(*.-^+'^.--^) 

(59) 



and from the first formula (58), 
K(t, s + [zi;])fc^£, with jw^-^l > \v 



1 



du 



\u\=p<l 



dv- 



1 1 



|„|=p-i>i 



u — n 1 



(60) 

We now apply this to the special locus L. Then picking p such that a < p <1 
and p~^ <\z\^ deduces from (59) and (60) the three kernels (49), 

dsTl) and dlol): 



\^{t - [z'\s)k/ 
K(t-[z-i],s + H)fe,^ 
K(t,s + [w\)k,i 



k+e T^(3) , 



kl ' 



(61) 



(-i)^'+^45(t.) 



In order to compute the polynomials (55), for given by (57), one needs 
to evaluate all the functions involved along the locus L; in particular, from 
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dSTl) and p4|, 



Z{t,s] 



(1 + a2)"("+i) / 0, 



and from (45), one finds for a < \z\ and \w\ < a ^, 

1 



Z(t,s+[z]) 



Z(t,s) L Z{t,s) 



e 



E°iiii.(«i+2Vi) 



(62) 



{1 + az) 



Z(t-[z-i],s) 



1 



Z(t,s) £ Z(t,s) 



e ^3 



= (1-- 



Z(t-[z-i],s + M) 



Z(t,s) 



(1-f) 



£ (l + ait;)"(l- f) 
(63) 



From the definitions of tlie Fredholm determinants (48), we liave from the 
expressions ((stI), ^ and (fsol) in terms K(t-[z-i];s) and K(t, s+[t(;]), 
that 



rfc(t,s)|^ = Z(t,s)/fi^)(0) = Z(t,s)//^^^(0) = Z{t,s) Hi"' (0,0). (64) 

Then, setting Hk{0) := i^P(O) = i?f^(0) = //f^(0,0), one finds for the 
orthonormal polynomials ( |55[ ) with regard to the inner-product ^ , ^= 



(2), 



r(3)/ 



( , )t,s , and for the Darboux sum (56): (remember notation (5)) 

Pk{z)=p^^\t,s;z)\^ 

fcZ(t- [z-i],s) det(l-K(t-[z-i],s))fc 



Z{t,s) y/det(l - K(t, s))fe det(l - K(t, s))k+i 



w 



,Z{t,s + [w]) det(l -K(t,s+ H))^ 



Z{t,s) ^det(l - K(t, s))k det(l - K(t, s)) 



fc+i 



(65) 



Y^P,{z)P,{w-') =J2pf\t,s;z) pf\t,s;w^')\^ 
j=0 j=0 

z\k Z{t- s + H) det(l - K(t - [z-%s + [w]))k 



w 



Z{t,s) 



det(l -K(t,s))fc+i 



establishing formulas (47), via (48), (^50^, m3\ and (64), with the double 
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(i) 

integral representation (49), (50) and (51) for Kf^ ^. Then using 

1 



1 1 



1 



1 {u — z) (v — w) 



1 



(z — w) 



(66) 



V — u {v — z) (u — w) V — u (v — z) {u — w)^ 



one notices that the kernels K^]{z ^) and Ku'^Uw) are rank one pertur 



(2), 



bations of the kernels Kj^f {0) and K'j^''j{0), whereas K^l{z '^,w) \s a rank 



(2), 



two perturbation of K^\{G), as given by the second equalities in (49), (50) 



and (51). The functions , appearing here, have each two different ex- 



u ^. The functions h^^ are 



pressions, one obtained from the other by u 

contour integrals, for which it is convenient, for later usj^ to include z in 
the contour, at the expense of introducing a residue; so, one has from the 



definitions of (53) 



-dv 



1 



+ 



ro 



^ ^ 2m{v - z) (-w)'=+Va(2n; v) (-z)'^'+Va(2n, z) 



1 



hf\w) 



(-z)'=+Va(2n,z) 
-dv Lpa{2n]v^^) 

hf\w) + {-wY+'M'^n,w). 



(67) 



+ i-wY+^iPai2n,w) 



- (i) . (i) 

Then the /i^ have the following form in terms of the , upon using the 
expansion l/{v - z) = v-^/{l -^ = 1 EoHV ^ 



dv 



ro 



2iri{v - z) (-w)'=+Va(2n; v) 



a=0 



dv ipa{2n,v 



27ri(-vY+^ 



v — w 



a=0 



-W 



establishing formulas (53). Upon using l/(w 



(68) 

'-'Y:'^{u/vr, the 



kernels ^^^^(O) and iqY(O) can also be expressed in terms of the g) s. 



■^it enables us to expand l/{v - z) = v-^/{l - f ) = ^ Er(f ''^ 
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oo 



yielding: 

<(0) = E^aUl = 42(0) (69) 
and therefore K^^^~^{0) = i^'^^^(O). This ends the proof of Lemma 



4.1 



It is convenient to introduce a new function $(2) := (2), depending on 
m/n and a, namely 

n/2 



el*W:=(-zr((l + az)(l-^))"' . (70) 



(i) 

From the decomposition (|49|) and (|50|) of the kernels Kj^ ^, it follows from 



(10) and (48) that, in particular, 



'2m+l 

where, by ([67|, 

£>2m+l 



(71) 



£>2m+l 



^)2"^+Va(2n; z) 

:5(i)(^-i) + r(i)(z-i)"~"' 



g-n4>(2;) 



a — Z 



(72) 



i?(^)(z):= Qf/.f(.) 

£>2m+l 



Qfhf\z) + (-z)2™+Va(2n; z) QS+.(- 

^>2m+l 1 

:5(2)(z) + r(2)(^)-''-^ 



-n<I>(2) ' 

with Q^'^ a function of £, with support [2m + 1, 00), defined by 

Q« = [(1 - X2„+ii^«(0)TX2™+i)-iX2„+i<7«](£), (73) 
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and 

e>2m+l e>2m+l 

^(1) (74) 
e>i ^ ' i>i 



and with (remembering from (|68|)) 

oo oo 



a=0 a=0 



So we have the following expressions for the orthonormal polynomials 
P2m+i and P2m+i on the circle, as in (47) : 



^--'-'">^\/isi--"'""<'+"'""(^-^'"<^ 



(76) 



In particular, the fact that the norm equals 1 implies an identity, which will 
be used later; namely by (39) and (25), find: 



^2m+i(0) / dz (l+az)'^ / a\n+i 

'1 [1 +az) 



1 — ^ P2m+1, P2m+1 ^ 
_ g2m+l(0) f _dz_ 
~ ^2m+2(0) /ro,, 2T,iz (1 _ 

(l-i?W(z-i)) {l-R^^\'/^ ^^^^ 

-H'2m+l(0) / 



(l-i?W(z-i)) (l-i?(^)(.)). 



i^2m+2(0) 7ro^ 27riz 
We end the section with some Fredholm determinant identities: 

Lemma 4.2 fl23^ ) Given a trace-class operator A and a rank one operator 
a®h, the following holds for an arbitrary constant c: 

det(/ - A + c a®h) = {I - c) det(/ - A) + cdet(/ - A + a®h). (78) 
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Lemma 4.3 ( \23^ ) Given two continuous functions F{z) andG{w), vectors 
a^, depending continuously on z and w, and a trace class kernel K, one 
has that 



detl I-K + 



+ 



F(z 



dz 



7ri 



Giw[ 



7r2 



dw \ 
w I 



det(/ - K) (79) 



5 Representations of the kernel IK^^^ for the 
domino outlier paths 



The purpose of this section is to deduce Theorem 1.1 from Proposition 3.1 
using the Fredholm determinant representation, either of the Darboux sum, 



or of the bi-orthogonal polynomials themselves, as in (47). They will lead 



to the two very different expressions for the kernel, announced in Theorem 

The following formula, based on contour integration, will be useful in 
the sequel: 



1 



{2m) 



dz 



ro., 



dw 



F{z,w) 



w — z 



{2-Kif 



dw 



ro,, 



dz 



ro,a,- 



F{z,w) 



w — z 



2ni 



(80) 



Fiz,z). 



5.1 Representation as a perturbation of the Krawtchouk ker- 
nel 

In this section we establish formula Q of Theorem |1.1[ To do so, we take 
formula (30) for the kernel ^n,m and the third formula of Darboux type 
in (47) as starting points. Recall some of the definitions given in ([4]); the 



definitions of ax,s and by^,^ will turn out to be the same as the one in (|4l). 



37 



Here we have some additional ones: 



Gy^r{w) 



2m 

{-i)y 



l+az\ 2 



1 



l + aw\ 2 



, Gy{w) = Gy^n/2{w), 



dz 



7ri 



dw 



GyAw)hf\w), by{e) = by^^/2{e). 



(81) 
(82) 



w 



We now extend the functions gf\k) of (52) as fohows: 

dz 



gf\k;s) 



gf\k-r) 



ro.a 



2-Ki 



{-zfipa{2k-z) 



1 + az 
1 - 



dz i-z)-^-"^ /!_+ az 
To,a 'Pa{2k,z) V 1 



a I ' 

z 



, g'l\k;0)=gl,'\k), 
gf\k;0)=gf\k), 



and gf\k,j), gf\k,j) are a smaU variation of g^^^(/i:,j), gf\k,j). 



(83) 



Proposition 5.1 The dual kernel has the expression (85). It is a perturba- 
tion of the one-Aztec diamond kernel ^ and involves the resolvent of the 



kernel (49 



12 



-lY-y^ZU'^T,x-2s,y) 
: -ts<r{-lT~''i^2{s-r) {x, v) + S(2r, x; 2s, y) 
+ {{l-Kill^,{0))-'a.y4k),b^,Ak)) 



(85) 



fi{2m+l,...) 



Compared to g^^', the expression g^^' has an extra term z ~ a in the numerator. Also 
the expression p^^' has an extra term w — a in the denominator; namely 



r(2) 



dz . ._^_2 z - — d 

— (-^) — -, — ^■ 



(84) 



'(/(fe),<7(fc))^2(2. 



m + l,...) Z^2m + 1 



/(fc)ff(fe). 
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with 



m 



S{2r,x;2s,y) 



/ 1/ 



'7r3 ^ •'7r2 ^ 



\27riy 

00 



7r3 •'7r2 



W — Z 



n n 



dzi — ^^v-^-. /l±a.\ 2/l±^ 



/ V y-y+m+b I 2 ' * 2 7 ^-a- 



(2) / n n 

+m+b I 2 ' 2 



O 1 O 



b=0 

and mi/i 



(86) 



(-1) 



(27ri)^ 



7=1 •^7ri 



Krii) 



■1) 



1+aui 



(27ri) 



^ I 111 — 7/1 in. 



'7=2 



u — w (pa{n; w)wy'^^- 



2 ' 2 



■^y-m+e V 2 ' 2 



6>0 



(87) 



Proof: Throughout the proof the fohowing radii, in between a and a~^, 
will be used: 

a < < r2 < S2 < si < ri < S3 < a^^. 



The contours Fo^a and TQ^a,z in the kernel (30) will be deformed to circles 



7^2 and 7^ respectively, with w z interchanged, 

{-ir-yKn,m{-y,-x) 



(27ri)2 



7r2 



It; 



dz 



z-^ 



2m 



p„ [Z] 

^■-1 fc=0 



(89) 



^^Note that the definition of S(2r,x;2s,y) is compatible with the one in Q. Also it 
appears as a part of the kernel ( 102 1 . 
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Using ( 47 ) , MSj) and (pTl) combined for the Darboux sum, the formulas ( 24 ) 



for p^{w) and Paiz), and the formula (46) for ipa{n; z), and Lemma 



4.2 



for 



the rank 2 perturbation, one finds for the integrand of kernel ( 89 ) 



{-if-y p^{w) pI:{z) z--^ 

{2Triy w z wy~" 
i-ir-y ipa{n;z) z-+'^ 



2m 



k=0 



det(l -Kj^l{z-\w] 



£2(2m,...) 



(27ri)2 ipa{n; w) u;2/+'"+i z-w det - Ki^}(0)] 

\ '"'^^ V<?2(2m+1,...) 

„ , , ,det(l- Ki^){0) + {z-w)h^\z-^)(-^hf\w))) 
F^{z)Gy{w) V k,n J \ J k \ n w e \ " J i2(^2m+i,...) 



^2(2m+l,...) 



F^{z)Gy{w) 



1 



z — w 



1 + 



det (l-i^g(0) + 4^)(z-i)(-l/.fH)) 



^2(2m+l,...) 



det(l-K^^^^(0) 



'£2(2m+l,...) 

Then one performs the if-contour integration over and the z-integration 



over ; from ( 80 ) , one has that for any a < rs < r2 < ri , 
dz Fx{z)Gy{w) 



dw 

dz 



z z — w 



dw 



^ Z z — w 



^ ^ ^ + 2m i — F^{w)Gy{w) 



with 



27Ti i —F^{w)Gy{w) 



J"/ w 



-lY-y r dw 



which is used in = in (90) below. Remembering the definitions (82) of a^ik) 



and by{i), this yields (using the radii (88)), thanks to Lemma 4.3 used in 



40 



equality = in (90), 



dw f dz Fx{z)Gy{w)w 



7r-2 



+ 



w 

dw 
w 

dw 
w 



7ri 



z — w 



dw 



7r2 



w 



Gy{w)w 



dz 



FJz) 



7ri 



dz 



Fx{z)Gy{w)w- 



det (l - i^g(O) + /,«(z-i)(-l/,f 



P{2m+1,...) 



dz Fx{z)Gy{w)w 



det (1 - kW(0) - a. 5j,)^2(2^+i_...) 



1 + 



z — w 



det(l -K(i)(0)) 



e^{2m+l,...) 



dz f dw Fx{z)Gy{w)w 



7r3 



7r2 



Z — 



x,y 



1 + 



det(l -i^(i)(0) -a^. (8)6j/) 



£2(2m+l,...) 



det(l-i^a)(0)),2(2„+i,...) 
Then the dual kernel ]K„^m(— 2/) —2;) is given by 



(90) 



7rc 



z 



dw Fx(z)Gy{w)w 



7r2 



w 



w — z 



+ 1 



det(l-K(i)(0)-a^(»6y)^ 



'(2m+l,...) 



det(l-KW(0)),2(2^^.i,..,) 



(91) 

thus establishing formula (85) for r = s = n/2. 

The second formula (86) for r = s = n/2 follows immediately from 
the definitions (81) of Fx{z) and Gy{w), while the third formula is proved 
exactly as in (68). It remains to prove the formulas (87) for ax{k) and by{i). 
Indeed, formula (80) permits us to interchange the integration variables in 
the expression below, so that \z/v < 1, enabling us to expand = 

Ylb>o(v^^- Thus, using the radii (88), and using the explicit expressions 



41 



(81) and (53) for and h!i \ one finds for ax{k) as in (87), but for s = n/2, 



ax{k) 



(-1) 



x+k 



(27ri)2 
(27ri)2 



z (Pa{2n; v)v'' v — z 
ipa{n; z)z^~^"^ 1 



dv 



(Pa{2n]v)v 



k+2 



1 



7ri 



^x+m-k-2 ^ ^ 

27ri V^aC"-;^) 



m+fc V 2 



, 6>0 



and, similarly for by{i), for r = n/2, 

/■ dw r du ipa{2n;u)u^+^ 1 



(27ri)2 



7r2 



y+m+2 Y 



(92) 



w — a 



, b>0 



This proves Proposition 5.1 for r = s = n/2. 

We now proceed to compute the extended kernel for general r and s, 
using the recipe (29) in Proposition 3.1 applied to the non-extended kernel 
([91]). Thus, 

(-l)^-^K-;,(2r,x;2.,y) 
= -l,<,(-l)"-^V2s-2r(a;,y) 

+ (-l)^-^V'n-2r(x, •) * K^^^(n, •; n, o) * V'2s-n(°, y) 

= -ls<r{-ir-'^^2s-2rix,y) + (-l)"-^' 



dz 



dw 



7r2 



W — Z 



{^n-2Ax, U) * G-u{vj){-ir) {F^,{z){-iy * i^2s-n{v, y)) 

+ (-1)--^ ((1 - i^«+i(0))-'a-.(A;)(-ir \ ^2s-n{v,y), 

'lpn-2rix,u) * &-« (^) (-l)")^2(2„+i,...) 
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= -l.<rV'2.-2r)(a^, y){-lT~^ + S{2r, x; 2s, y) 

We now show the last identity in ( |93[ ). The useful tool here will be the 
following formula for the *-product (with respect to x) holds: 



Z * 



J 2mw ^ ' J2Tnw^\wJ ^ ' ^ 
from which it follows that (see Ml)) 



(94) 



Z " % tp2s-niv,y) = Z' 



y(\±az\ 2 



1-^ 



(95) 

We thus need to compute the following *-product, with regard to u, also 
using formula (|95|), 



(-l)^F„,,,(z) = F„4z)(-1)^ * ^2s-n{v,y) 
_ y^a{n; z) 



2TTi 

V>a{n;z) 
2iTi 



-y+m+l j l+az 



* '^2s-n{v,y) 
s-n/2 



(-l)^F_,(z)(i±f 



-n/2 



(96) 



from which one deduces the following *-product, also in v, remembering the 
definitions (53) of the h, , and using the definitions (81) and (82), 



dz _ 



7ri 



dz 



h^^\z-^)F^,{z){-lY*^2s-n{v,y) 

F_,,,(z)(-l)^4^)(z-i) 



(97) 



Similarly, one checks, using (81) and formula (95), 

i-iyC^riw) := lPn-2r{x,u) * G_4l«)(-1)" 

1 



-0„_2r(x,n) * 



1 w' 



.x—m—1 



27ri (pa{n; w) 

n 

1 + aw\ 2 ~' 



2-ITi ifain, w) \ 1 



1 - ^ / 



(98) 
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from which one deduces, using (IsTl) and (82) 



dw 



G. 



7r2 



The last identity in (93) now fohows from (96), (99) and (86). The for- 
mulas (86) and (87), involving the g^\^, •) are obtained exactly as in the 
computation for r = s = n/2. 

The identification of the first line of the equality ( 85 ) with the one- Aztec 
diamond kernel (|2|) is immediate, after setting z i— )• —z~^ in the integral 



V'2(s-r)- This concludes the proof of Proposition 5.1 



5.2 Representation of K^""^ ^ double integral 



Here we show formula (10) for K^^^ in Theorem 1.1 and also, along the way, 
an alternative representation (102), which will be useful in taking saddle 
point limits. The starting point here is again formula (30) for the kernel 
IC^^rra as in previous section; but instead we will now use the Christoffel- 
Darboux formula and the two first formulas of (47). The methods and the 
notation in this section are very close to those used in [2]. Remembering the 
definition (iTOl) of ^(z), define the functions 



Ei{z,w) := 62 

n 
62 



(l-5«(z-l)) (l-5(2)(«;)) 



W) 



E2{z,w) := et('fW+*(-))(z,; - a) {l - S^'\z-')) T^^\ 

E,{z,w) := e-t(*(-)+*(-))— r(i)(z-i) (l- S'^'Hw)) 

z — a \ / 

Ei{z,w) := 



and 



C{u; x) := Ci{u; x) + 2C2(n; x) 



27ri 

2 

~ 27ri 



(100) 



-^(l-i?«(z-i))(l-i?(^)(.))(i±|)" 



(101) 
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Remembering the form of the kernel ]K„^m(— x, — y), obtained in (30), one 
has the following Proposition: 

Proposition 5.2 The dual kernel K'^^(2r, —x; 2s, —y) takes on two differ- 
ent forms 

(i) the form (10), as stated in Theorem \l.l\ 

(a) or another form, useful for taking saddle point limits (see section 7.2): 
SSfM^'^-(2r, -x; 2., -y) = C{s-r;x- y) 



ls<r(-l) 
/ 



X-y H2m + 2iO) I , 



dz 



+ 



0. (2vri)^ 

n 

1-J 



dw 



\ 



YXEi{z,w) 



z — w 



V 



l-\-aw \ 
l_iL ) 
w / 



-w 



\x+l 



l±az \ 2 



(102) 



Proof: As a first step, we prove the case r = s = n/2. The kernel M.n^m{x, y) 
in ( 30 ) can be expressed as 

1 f dz „, ^ r dw .w-y-"" 1 



(103) 

X ((f)'™+'P2m+l(^)P2m+l(u;"') -P2m+l(^"')P2m+l(u;)) , 



using the Christoffel-Darboux formulsp] for bi-orthonormal polynomials on 
the circle, (for M = 2m + 1) 

E P.(.-')n(») = '-■'Pm(.K'PmC.-')-Pm(^-)Fm(«.)^ 



A:=0 



^^This can be shown by generalizing an argument of B. Simon in |27) . first proof of 
Theorem 2.2.7. 
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One then uses the representation (47) of the polynomials in terms of the 
nj^^^s, yielding: 

dz 



dw I {l + az) (1- f) 



n/2 



W 



-y—m—l -j^ 



ro,. /ro,.,. 27ri Ul + aw) (l - £ 



z — w 



H. 



(1) 

2m+l 



-ff2m+l(0)-ff2m+2(0) 



ro,, 



/■ / (l + a^;) (1- ^) 



n/2 



1 - ^ \ W 



Z 2:+m+l 2; — IL) 



-ff2m+l(0)i^2m+2(0) 

(105) 

Using (80), one finds, using (71) and (72) for the Fredholm determinants 
and the function ^{z) as in (70), and setting x i— )• —x, y i— )• —y, 

-ff2m+2(0) 



H2m+l{0) 



^n,m{ 27, y) 



V27ri/ /ro 



X (l-i?(i)(z-i))(l-i?(2)(^)) 
dz z^"^^^ /l-r;\ S{$(u,)-<i.(^)) 



w — z w y 



1 



e2 



X (1 -i?«(w-i))(l 



+ 



27ri 



dz 



ro,a 



^(l-i?«(z-i))(l-i?(^)(z)). 



(106) 
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Using the identity (77) applied to the last term of Kn^mi—x, —y), the dual 



kernel is then given by: 

-y) = — " ^n,mK-X, -y)) 



-f^2m+l(0) 



J Vc\n. J Fn 



2vri Jr^ ^ 



(107) 

Multiplying the kernel (107) above by (—1)^ ^, interchanging z ■h^ w m. the 



second double integration and then combining the two first integrals, one 
finds that the kernel 'Kn,m{—x, —y) consists of two parts, a double integral 



3C(x,7/)and a single integral Ci(0;a; — y), defined in (101), 



(-l)"-" '^'""''i°i ]K„,^(-^, -y) = y) + Ci(0; X - y), (108) 



where (recall -R*-*^ is given by (72)) 



dz 



y-i ^ 1-1 



X et(*W-*(-))(l - i?(i)(z-i))(l - R^^\w)) 

(109) 



and, confirming (101), 



Ci(0,x-y) 



2m 



dz 



To. 



(l-i?W(z-i))(l-i?(2)(z)). 



Notice this establishes formula (10) of Theorem |1.1| for r = s = n/2 
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Using (100), the double integral-part % of the kernel (108) reads, upon 
multiplying out (1 - i?(i)(z~i))(l - R^'^\w)), 

%{x,y) 



{2m) 



dw 



1 



a — w 



62 ' 



1 _ - r«(z-i) 1 ( 1 - 



a — z 



-n<I>(-!ii) 



T(2)(ti;) 



(27ri)2 



i) dw- 



Z — W 



, (-z)^ i-f 



+ 



i-zY {-wY+^ 1 - ^ 

(110) 

The part of the double integral, involving £4(11!, z), is not in a usable form, 
in view of the saddle point method and the topology of the contours; the 
integrations have to be interchanged, at the expense of a new residue term. 
So, using (80), the double integral %{x, y) involving E^lw, z) becomes, after 
a further interchange z ^ w: 



(27rz) 



dz 



0,a 



dw ({-z)y-^ {-wf - z I - 



+ 



T,a..Z-W\{-wY {-zY+^l-tJw 



dz 



—^T^^\z-^)T^^\z). 



/ro,. {-zY-y' 

Defining L{x,y) as the double integral part of %{x,y), namely 



(111) 



£(x,y) := 



(27ri)^ 



dz 



dw- 



TXE^{z,w) ({-w)y-^ , {-z)y i-f 



0,a •'J-0,a,z 



Z — W 



+ 



i-zY i-wY+'i-^ 

(112) 

and remembering the definition (101) of 6*2(0, x), formulas (110) and (111) 
show that 

%{x,y) = L{x,y) + 2C2{0;x - y). 
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Combining this formula with (108), together with (101), yields 
-ir-y?^^^^^Kn,m{-x, -y) = %{x, y) + Ci(0; x - y) 



2m+l 



(0) 



= y) + 2C72(0; x - y) + Ci(0; x - y) 
= ^ix,y) + C(0;x - y), 

(113) 

thus establishing formula (102) of Proposition 5.2, with the Ei and C(0;x) 
as in (100) and (101), in the case r = s = n/2. 



The proof of the extended case follows the recipe (29) in Proposition 3.1 



applied to the non-extended kernel ]K„^m(x,y), namely: 



{-ir-y^P^,K'^^U2r,x;2s,y) 



1 



s<r\ 



-IT ^i^2{s-r){x,y) 



H' 



2m+2 



(0) 



2m+l 



(0) 



+ Vn-2r(x, •) * (-l)"-^ ^'"^+'i°i lKn,^(- ;o)*V'2.-n(o ,y). 

-n2m+l(Uj 

Notice that the *-product only acts on the exponents of w and z in the 
expressions (109), (112) and (101) for %, L and Ci(0;x). To do so, one first 
removes the —signs in —z and —w and one changes x i— )■ —x and y i— )• —y 
and ones use the identities (95) below: 

-v-l 



'lpn-2r{x,u) \ 



W 



+ 



1 



Z 



-u+1 I 



T % ll^2s-niv,y) 



w 



+ 



1+aw \ 2 

1-^ 



1+aw ^ 2 



l+az A 2 
1_» 



and 



V'„-2r(2;,u) *„( — )*„ i'2s-n{v,y) = Z' 



x-y 



1 + az 

1 - 



(114) 



One now applies these formulas to the two expressions (108) and (113) for 
the non-extended kernel 'Kn^m{x^y): 



(-ir 



i{x,y) 



-^ -^2m+2(0) j 

^2,n+l(0) (115) 

%{-x, -y) + C7i(0; -x + y) = £(-x, -y) + C7(0; -x + y). 
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where % and L are defined in (109) and (112). From (114), the *- 
multiplication amounts to inserting some elementary fractions. Thus the 
irst repr( 



first representation (115) leads to 



+ 



+ 



H2m+l{0) 
-ls<r{-lT~^^p2(s~r)ix,y) 
5x^y f dz 

2m 

' ' ) 



{2r,x;2s,y) 

H2m+2{0) 



H2m+l{0) 



\27riJ 
( 



(p dz (p 



dw 



W 



(1 - i?(i)(z-i))(l - i?(2)(^))et(*W-*(-)) 



i-w) 



1+02^ 2 

1-2 



V 



i-z)- 



n 

1+aw \ 2 
1_A 



+ 



{-z)-y 



l+aw \ 2 
1_A 



{—w) ^ 



+ 1 



/ 1+az 



J 



(116) 

which is expression (i) for the extended kernel (formula (10)), upon using 
the explicit expression (fral for e'^*^^)/^^ Expr ession (ii) (formula (Il02|) is an 
immediate consequence of the second representation ( 115 ) for 'Kn^mi—x, —y), 
with £ as in ( 112 ). This ends the proof of Proposition 5.2 , and also the proof 
of Theorem II. 1[ ■ 



6 Scaling limits of various functions 
6.1 Limit of the functions (7^*^ 



Remember the quantities fo,^, p, 6 defined in (12) in terms of the weight 
< a < 1 on vertical dominoes. Also the reader is reminded of the scaling 
(jTsJ): 

n = 2t, m= r + apt^/^ 

a + a-i ^ (117) 

X = 2a^eTt^/^ + s = t + (1 + a2)^rt2/3. 

In the expression (j7|, the extended kernel involves discrete variables 

k and £ and involves integration variables u and v in i^«(z-i), which we 
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rescale as follows 
At 



k 



a + a 



At 



U = Vq + 

z = vo + 



U 



A{2ty/^' 
c 



a + 

V = Vo + 

w = vq + 



:j + Ap(2t)i/3 + 1 



V 



A{2ty/3 



(118) 



The va riabl es z and with the rescaling above, will appear much later in 
Section 6.3 Also note that k,i > 2m + 1 implies k, X > a := 2'^^^a. 

Define the function F{v) and G{v), together with their Taylor series 
around the saddle point v = vq := — and the Taylor series of log{—v) 
about Vq, 



F{v) := log(l + av) + log(l 



- H 



:jlog{-v) 



G{v) 



-- (1 + a')e log - 2a^e log(-?;) 

V 

= G{vo)-A\v-vof + 0{v-vof, 
log(-w) = log(-t;o) + :^{v-vo) + (D{v - ^;o)^ 



(119) 



with 



(120) 



G{vo) = il-a^)9log{-vo). 
Given parameters r and A, set for future use: 

MX(t) := 2tF{vo) + rG(7;o)(2t)2/3 + Xplog{-vo){2t)^/^ 

= 2tF{vo) + ((1 - a2)0T(2t)2/3 + pA(2t)i/3) log(-z;o) (121) 



M,(t) := Mlit) 



r=0 



Note that, upon introducing the scaling ( 118) for u and v, the function ^{z) 

as defined in (70), has the f ollow ing Taylor series about vq, which then can 

V 



be scaled using the scaling ( 118 ) for w — 



n^{v) = 2tF{v) + ap{2t)^^^log{-v) 

= Mait) + -aV + Oit-^^^), 



(122) 



51 



where F and M5- are define d in (119) and (121). Combining this with 
log{—v) and, using the scahng nm and mn for £-2m = {X-d)p{2ty/^+l, 
and V — vq, yields: 

n<!>{v) + {£- 2m) log{-v) = 2tF{v) + Xp{2t)^/^ log(-u) + log(-w) 

= Mx{t) + log(-i;o) + -XV + 0(t-i/3). 

(123) 

We s hall need the following expressions in terms of the scaling (117) and 
(118), together with the summation variable b = (3p{2t)^/^ : 



k + b 



At 



a + a' 
2t 



:^ + {K + p)p{2t)^/^ + l 



-x + m + b= - 2a20rt2/3 + (a - ^ + 2^/'^^)pt^/^ 

a + a ^ 

2f 

x-m + £= + 2a20rt2/3 + + ^ + 2^/^X)pt^/^ + 1. 



(124) 



To do the asymptotics, the following Lemma from Borodin-Ferrari [10] comes 
in useful: 



Lemma 6.1 Given a contour integral of the form 



* 2m 



(125) 



7{o,pj,...,p„} 



where 7{o,pi,...,p„} cl closed contour, having as poles the points 0,pi, . . . ,Pn 
only. Assume that ^fo{z) has a real saddle point vq and a steepest descent 
path^' about Vq; assume i/io^ 7{o,pi,...,p„} can be deformed toj', without pass- 
ing through the points {0,pi, Let fo{z) experience a double critical 
point at Vq, and fi{z) a single critical point: 

foM = foivo) = 0, f^'ivo) > 0, fiivo) = 0. 



Then, remembering the definition (I4) of 

lim^_(fACo)i/3eT[*/o(-o)+t^/='/iK)+t^''^/2(.o)+/3(.o)]j±=sgn(„o)Ai(±-i)(K2), 

(126) 



where the convergence is uniform for ki, K2 in a bounded set, and with 



Kq 



2k, 



2/3 ' 



K'2 



.1/3 ■ 



(127) 
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For any choice of L > 0, 



uniformly for K2 > —L and ki in a bounded set; cq and ci are positive 
constants independent oft. 



(i) 

The next Lemma deals with the scahng hmit of the 5^ s, given by (52) 

I I I I I 1^ 

and (83) and the 5^ s, given by (84); also remember expression (121) for 
Mx{t). 



Lemma 6.2 (Limits of 5^*^) Given the scaling (118) for I, together with 
b = /3p{2ty/^, the following limits hol^^ : 



(129) 



lim (2t)V3^^e-M,W a)(2i) = _Ai(A) 
t-5-oo vq — a 

lim {2tY/^A{vo - a)vle^^'^''^gf\2t) = -Ai(A), 



t— >-oo 



where the limit is uniform for X in a bounded set. Also, for I as above and 
for any choice of L > 0, 



(2t)i/3eTM.(t)^(2)(2i) 



< CqC 



(130) 



uniformly for A > —L, with positive constants cq and ci independent of t. 



Also, using the scaling (in) and (118), one obtains 

lim tV3^,-*C.V3,(*/2)^(i) ^^.(_.)(^_^^2V3^) 

t-s>oo Vq — a x-r/(i.-ru 

lim t'/'Avo{vo - a)e'''U^-^'',^''^\^^l it-t- s) = Ai^^^a - ^ + 2^/3/3), 

(131) 

where the limit is uniform for ^, r, /? in bounded sets . Also, for any choice 
ofL>0, 



< coe 



(132) 
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uniformly for — ^ + 2^/^/3 > —L, with cq and ci positive t-independent con- 



stants. Furthermore, with the scaling (111) and (118), 



lim tV3^g-^C+,+.V3,(*/2)~(i) ^ _ ^ _AiM(-a + ^ + 2V3a) 

mn tV3^^2^^::,,,,V3.(*/2)^2)^_^^(^. ^ _ ^) ^ _Ai(-)(-a + e + 2^3 A), 



(133) 

where the limit is uniform for X in a bounded set and for any choice of 
L>0, 



uniformly for ^+2^/^\ > —L, for cq and ci positive t-independent constants. 
Proof: Indeed, in view of a saddle point argument, consider the function 



^F{x + iy), for F{v) as in (119); it tends to — oo near {x,y) = (a,0) and 



{x,y) = (— l/a,0) and tends to oo near {x,y) = (0,0) and at oo on the 
Riemann sphere; near its saddle point v = vq := < 0, the function 

?fiF{x + iy) looks like 3f?(x + iy)"^, with the three separatrices through v = vq 
in the upper-half plane y > corresponding to the three separatrices of the 
3^(x + iy)^ profile in y > 0; they are indicated by dotted lines in Figure 18. 
So, the levels profile consists of (topological) circles about 0, a, —a~^ and 
oo, separated by the separatrices. A path C of steep descent departing from 
the saddle point will be given by, say, extending a small straight line making 
an angle tt/Q < 6 < 7r/2, as indicated in figure 18, until it hits a level curve 
in that quadrangle, which one then follows until it hits the x-axis. Then 
one takes the mirror-image with respect to the x-axis. The contour thus 
obtained winds around 0, but not around a. So, for the function qi^\ one is 



in the condition of Lemma 16. 1[ 

The function g^^^ , involving integration about a contour Tq (about and 



not about a), takes on the following form, using the first equality of (123), 



J^^(2t) = / Z^^(l_«)e2**(^)+(^-2'")log(-^') = Zf^g2t/o(i>)+(2t)i/3/2(„)+/3(„)^ 

with 



foiv) = Fiv), f2{v) = \p\og{-v), f-siv) = logi-v) + log(l - f ). 

with 

/oM = /o (^o) = 0, y[^'ivo) = A^>0, f^ivo) = -AX<0, 
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satisfying the requirements of Lemma 6.1 

Then, by Lemma |6.1[ one has that for any L > 0, keeping in mind the 



definition (121) of Ma (t) 



(2t)V35W(2t)e 



-2tfo(vo)~{2ty/^f2{vo) 



-ci A 



(135) 

uniformly for A > —L, with positive constants cq, ci, independent of t. 



Then from (123) and using the usual saddle argument, it follows that for 



A E K, a bounded set, 



(2t) 



l/3„-MA(i)„(l) 



a-VQ 
A 



dV 
y 2-Ki 



with 0{t~^/'^) just depending on the set K. The uniformity follows from 
Lemma [6.1[ Indeed, in the limit t — >• oo, the only contribution of the integral 
comes from a small neighborhood of the path C of steepest descent, near the 
saddle point vq in ^-coordinates; this yields two rays emanating from in 
V-coordinates with an angle 7r/6 < 9 < Tr/2. This proves the uniform limit 
^ for 9f)(2t). 

The limit (fl^gl) for 



gf\2t) 



-dv g-^^W-C^-Sm) \og{-v) 



2m 



is similar, except that, since n^{v) + — 2m) log(— f ) appears with a minus 
sign in the exponential, one looks for a path of steepest ascent. Such a path 
consists of a loop about and a, passing through the saddle: take a line in 
the direction 7r/2 < < 2tt/3, until it hits a level curve, then follow this 
level around a; then, take the mirror image with respect to the x-axis. So, 
this leads to a steepest ascent contour about and a, as depicted in Figure 
18. Then, by Lemma 6.1 or more directly upon taking the limit, when 
t — 7- oo, the only contribution comes again from a very small neighborhood 
of Vq. One also has a similar uniform estimate by the arguments of Lemma 
6.1, for any L > 0, A > —L, 



(201/3^(2) (2t)e2*/oK)+A(2t)V«/2(.o) ^ (2t)l/3^(2)(2t)eMA(t) 
with t- independent positive constants Cj. 



< coe 



(136) 
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Figure 18. Plot of level curves^® for ^F[x + iy), for F[v) as in (119) 
with a = 1/2 and saddle at vq = = —1/3. The six dotted curves 

(departing from vq) are the separatrices of the level profile of + iy) 

and correspond in the lower picture to the six separatrix curves of ^{x + iy)^ 
departing from the origin. 



^Near very dark regions the function tends to oo and near clear regions to ~oo 
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To prove ( 131 ), we invok e Lem ma 6.1 applied to the following functions, 
with the scaling (117) and (118), together with b = I3p{2t)^/'^. So we must 
prove that 

s-t) = -I+{t) and - s) = -/"(t) 



with of the form (125). Using the explicit expression (83) of g^^\t;s — t), 
using expression (|122D for $(z), formula ploj ) for G{z), and the scaling ([124]) 
for —X + m + 6, one checks : 



l+az 



^$(2) + (b_a.)l0g(-^) + (s-t)l0g^+l0g(l-|) 



ro,a 



ro,a 



271-i 



2^' 

,(2) 



and similarly for g^^x+m+bi^'^ s — t). So, for both cases: 

foiz) = F{z), h{z) = TG{z), /2(^) = (^-e + 2i/3/3)plog(-z) 
log(l-f), for/+ 
^ log((l-f)z2), forlr ' 

for which one checks from ( |119 ) and (12) that 

Ko = A^ Ki = -T, K2 = a -^ + 2^/^13, 



h{z) 



and consequently, from (121), 



log^ 
\og{vQ-a)vo 



This and sgn vq = —1 leads to the limit (131), together with the bound 
p2l). 



To see (133), one shows that 



with /.^ as in Lemma 



(t; s-t)= /+ and £l^_,,{t, s - t) = 



6.1 



Using (124), one checks, for both cases, 



fo{z) = F{z), h{z) = -tG{z), f2{z) = {-a + C + 2^/^ \)plog{-z) 
h{z) 



0, for /+ 



log z , for 



57 



for which one checks that 

Ko = £', Ki = r, K2 = -a + i + 2^'^\. 
and consequently 



This leads to the limit ( 133 ) and the bound ( 134 ) , ending the proof of Lemma 



6.2 Limit of a_y^s{k), h^x,r{^)i K^u and S{2r,x;2s,y) 



In the usual scaling (117) and (118) one replaces x by two variables x and 



y, and also s by r and s, leading to new variables ^i, ^2 and ti, T2 ; so the 
scaling is now: 



scaling (117) and (118), together with 

X = 2a'6nt^/-' + Cipt'/^ y = la^Qr^t^l'' + i^pt^l'' (137) 
r = t + (1 + a^)eTxt^l^, s = t + (1 + a^)eT2t^'^. 

The reader is reminded of the definitions (14) for A^^\k) and 
47'-^^)(ni,n2). 

Lemma 6.3 (Scaling limit of a^y^s{k), b-x,rW and Kj^^{0)) Given the scal- 
ing ofx,y and k,i,r,s as in (131), one obtains for a-y,s{k) and b-x,rii), 



A/^^'^ (t/2) 

hm i-iripvot'/'') e e2-+2V3.^/ ' a.y,s{k) = A-J^Jk) 

1 /q -Af"^l , ,,, (t/2) 



(138) 



with the limit uniform for ^i, 6.2, ti, T2 in a bounded set. Also for any L > 0, 

-ci(6+2l/3«;) 



1 /q M ^ , (t/2) , ^ 



,0 -M ^1 {t/2) 
-/3g _<t+5^+21/3a^ ' 



< coe 



< coe 



^ci(6+2V3a) 



(139) 
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uniformly for £^2+^^^^^ and £^i+2^/^\ > —L forcQ, ci positive t-independent 
constants. Finally, 

^lim (2t)i/3e*^-(*)i^(;)(0)e-^'^-(*) = ^lim (2t)i/3e*^-Wirg)(0)e-*^'^W 

f oo 



d/3 Ai(K + /3)Ai(A + /3) = ifAi(K, A). 

(140) 

where the convergence is uniform for n, A in bounded sets. One has that 



(2t)i/3eA^A(t)^(2)(o)e-A^4*) 



< coe 



-ci(k+A) 



(141) 



uniformly for k, A > —L, with cq and ci positive t-independent constants. 



Proof: In the limits below, the M's defined in (121) satisfy the following 
identity, useful later on: 

K-^+2^f^pit/^) + ^C+€+2V3 = M,+p{t). (142) 

Setting n = 2t, as usual , and using ([I29j), ( pT| ) and pl2| ), one checks, by 
substituting the scaling ( 124 ) for —x + m + b and k + b, with b = /?/3(2t)^/^, 
that 



lim pvot^^^ e"-<T+ii2+2^/^^ 

t—^oo 



M 



(i/2)v^ (1) 



6=0 



E 



(At) 



(2t)V3A(a-.o).o^e^-^W5f) 



/>oo 

/ Ai(-^2) (a - ^2 + 2^/3^) Ai(K + I3)d(3. 
Jo 



(143) 

The last limit holds by dominated convergence, using the estimates ( |132 ) 
and (130); the sum above does not exceed in absolute value the following 
sum 

(2t)-V3co Yl e-^i('^-«2+2'^'^)e-'=i('^+^) < c'^e-'^^-^'+'^l 



Also, from (133), 



hrn (pt;oti/3) e tln^+ki^; s - t) 



t—>-co 



(144) 
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Combining (143), (144) and (87), one concludes the first formula of (138). 
From (129), (131) and (142), it follows in a similar fashion that 



lim y 



6=0 

A 



(2t)i/3^^e-MA+,(*)<^W^(2t) 



a-vo 



f 

Jo 



Ai(^i)(a - a + 2^l^l3)ki{p + \)dp. 

(145) 



From (133), one obtains 



:i/3j*/2)^(l) 



(t;i-r) = Ai("i)(-f7 + a + 2i/3A). 

(146) 

Assembling (145) and (146), one deduces the second formula (138) from 
(87). Next one uses (130), (132), (134) and an argument analogous to the 
previous argument showing the estimate for dominated convergence; one 
deduces (139). Finally, a-y^s{k) and b-x^ri'd) tend to their respective limits 
uniformly on bounded sets by the "uniform convergence on bounded sets"- 
statement in Lemma l6.ll 

(i) 

In order to take the scaling limit of the kernels Kj,:, using the scaling 



(137), one sets b = p/3(2t)^/^, and uses the infinitesimal elements 



Ai 



/9(2t)l/3 



dX and 



Ab 



p(2t)i/3 



dp, 



(147) 



where A£ = 1 and Ab = 1; using these, one obtains from (54) and (129), 
nme^-^'=W4J)(0)e-^-^^(*)A^ 



t— >-oo 



hm ^ ^/,„ 
t^oo (2t)V3p 



^ Ab 

^ , (2t)V3^ 



A 



VQ-a 
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dX / d/3 Ai{K + ^)Ai(A + /?) = i^Ai(K, A)(iA. 

JO 



(148) 



Note that the sum in = can be viewed as an integral of a piecewise contin- 
uous function. Then Hmj_j.oo and inte grat ion c an be exchanged, by using 
dominated convergence. Indeed, using (135) and (136), the sum in =, before 



taking the hmit, denoted (say) by {k,X), can be bounded in absolute 
value by 



xi'\K,X)\<co 



1 



^ , (2t)i/3 



(149) 



for X, K > —L and with positive constants q independent of t, yielding ( 141 ) 



This justifies the existence of the uniform limit to the Airy kernel for k and 
A in a bounded set. Moreover, since 

(2t)i/3gA/.(t)^(2)(0)e-M4*) = (2t)i/3eM,(t)^(i)(Q)g-A./4*)^ 



one deduces from the above 



hm (2t)V3e*^A(*)^(2)(0)e-*^«(*) = hm (2t)V3eA^A(i)^a)(o)g-Af4*) = KAi{X, k 

(150) 



t—>-co 



This concludes the proof of lemma 6.3 



Lemma 6.4 (Scaling limit of S{2r, x; 2s, y)). One obtains for the scaling of 



X, y, r, s as in (131), 

hm pti/3(-z;o)^-"+^-^5(2r, x; 2., y) = K^^r''\a - 6, - 6), (151) 
with uniform convergence for ^1,(^2, ti,T2 in bounded sets. 



Proof: At first, one checks that the scaling (137) and the definition (121) of 
imply 



(y _ ^ + ^ _ s) log(-r;o) = ((1 - a^)e{n - T2)t"^ + /9(6 " 6)*'/') log( 

= M^7^(i)-M^7Hi) 

(152) 
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Set, as usual, the scaling (137), with n = 2t, in the third expression (86) for 
S{2r,x; 2s, y). However, we now set b = j3pt^/^ instead of b = I3p{2t)^^'^; i.e., 
replace 2^/^P by (3. Then, one checks, using the limits (^13lL 



lim ptV3,( V(*/2)-Mf;^(*/2))^(2r, x; 2s, y) 



lim > 

t-5>00 ^-^ 



Ab 



/3 -MZ,^_^,{t/2) (1) 



— e 



9-x2+m+b 



{t; s - t) 



/3 



roo 

= / Ai(--2)(^_^2 + /3)Ai(-i)(a-6 + /3)d/3 = A'i?'~"'^(cT-ei,^-6), 

where we have used the exponential estimates (132) to get dominated con- 
vergence in the above limit to the integral in the style of (149). The conver- 
gence is uniform on bounded sets by uniform convergence on bounded sets 
in Lemma 6.1 and dominated convergence in the integral, concluding the 



proof of Lemma 6.4 



Lemma 6.5 Upon using the scaling of ^13^, with t\—T2 > 0, the following 
limit holds for the function ip2{s-r){x ^ v) ^'^ 

(ei-S2)^ 



lim pty\-Vor-^+'-\-lT-y^2is^r){^,y) 
t—>-oo ^ ' 



n - T2) 



with uniform convergence for £^i,^2-,ti,T2 in bounded sets. 



Proof: Using the scaling (137) in (Q and (119) 

dz 



ro.a 



27ri2; 



(-^) 



x-y 



1 + az 

1 - 



To 



2TTi 



where we replaced ro,a by Lq, since s — r < 0, and with 

fo{z) = {T2-n)G{z), /i(z) = (6-6)pln(-z), f2{z) = -H-z). 

By ( 119| ), G{z) experiences a simple critical point at vq, which suggests a 
saddle point analysis about vq; therefore, set 



Z= V2(ri-r2) Aiz-VQ^^^^ 



62 



and picking a path of steepest descent, one finds that 



It = - , ( / —e'^ ^^(^1^^^ + 0(^-1/3 

-Avo^j2{n - T2) \J-ioo 



I ^47r(Ti - ra) 



with 0(t^i/^) uniform for ■^a, "T"!, ^"2 in bounded sets. The proof is com- 
pleted upon exhibiting a path of steep descent through vq for the function 
G{z), or equivalently G{z). 



G{z) 
0(1 + a2 



= log(2; + a 1) - log(2; - a) + f ° , ° ) log(-2;) 

\a ^ + ay 

and note that ^^(^(z) at {—a~^,0,a,oo) is (— oo, — c«, +oo, +oo). To con- 
struct a steepest ascent path, through vq, start with a small line segment 
through vq making a right angle with the real axis and join it to a level 
curve in the upper half plane which it meets and follow the level curve until 
it hits the x-axis to the right of a; then reflect this entire path in the x-axis 
to complete the path. The construction of such level curves is described in 



the proof of Lemma 6.2 This concludes the proof of Lemma 6.5 



6.3 Limit of the ingredients of the integral representation of 



Recall the scaling (117) and (118), in particular the one of z and w. Also 
remember the following expressions from ( 16 ) , ( 73 ) and also define the kernel 

Q{k) = [(1 - xaKAiX^r^aAi] (k) 
= [(l-X2™+iK«(0)TX2„+i)-iX2„+i5»](^), 
3Cf (k,A) := (2t)V3e-*^«Wi^(2)(o)eA'^^W 



scaling ||118| 



One shows the following Lemma: 
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Lemma 6.6 The following limits hold: 
lim (2t)i/3^^e-^«WQl'^ = lim {2tY'^A{vQ - a)vle^'-^''^Qf'^ = -Q{k) 

t—^oo Vq — a t—^oo 

(153) 



fOO 

lim {vo - a)voe^^^^'^h^\z~^) = - / e-^'^'^^Ai{K + /3)d/3, 

(154) 



lim {{vo - a)vo)-'e-^'-^%f\z) = - / e^'^'"'^^ Ai{f, + P)d/3, 
lim {vo - a)-ie-^^^Wr(i)(z-i) = -e'^^V^^ / Q(^)gK2V3^^^ ^ -e-2-?Q(C). 



lim (t;o - a)e^-Wr(2)(z) = -e'^^^'^^ / Q(K)e-'^2i/3c^^ ^ -e2-?Q(-C). 

(155) 

lim S^^Xz-'^) = ^(C), lim 5(2)(z) = ?>(-() (156) 

w/iere 

j>{() = - Q{n)dK e-^ '^^Ai{K + /3)dp, 



with uniform convergence for k and C in bounded sets. 



(i) 

Proof: Limits of Qj^ : The proof of (153) proceeds along similar 
lines as in |4j, namely first one notices that the operator norm of the 

(2) 

Airy kernel ||x5--f^Ai(A, k)x5-|| < 1> and since Xd-^t i'^j^)Xa tends point- 
wise to X5-f^Ai(A, k)xct (with exponential domination as in (149)), one also 

(2) 

has IXa^t {i^j '^)X^\\ < 1 for large enough t from which it follows that 
1(1 - XaX?\K, A)xa)"^|| < C uniformly for t large enough. Because of the 
estimate (135) involving g£{l){2t), one checks that 

satisfies 

\Qt{i^)\ < ce-'^ 
for any > and some constant c > 0. 
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We have, using (140), (129) and the above discussion, 

hm Qt{K) 



t^co vq — a 



t— >oo 



hm ( 1 — e 

t— 5>00 



= [(1 - XaKMXa)'\a^ (k) = -Q(k), 

uniformly k in a bounded set. Similarly 



-1 



A{2t) 



1/3 



VQ-a 



(2t) 



lim (2t)i/3A(t;o - a)t;ge*^«WQ 



t— >oo 



-Q{k). 



(157) 



(158) 



Limits of /i^*^ : With the scaling b = /9/3(2t)^/'^, using the infinitesimal 
element Aft = 1 as in (147), and using the scaling (118) for z and fc, one 
checks from ([53]), ([129]) and that 



lim (7;o - a)foe*^''(*) V( 



fe=0 



lim (uo - a)uo 1 



lim 



b=0 



p(2f)i/3; 

45,(2t)V3pd/3 

-2^/^C/3Ai(At + /3)d^, 



using dominated convergence from (130). Similarly, 

) 

^""^PM{k + /3)d/3. 



lim ^ ^ — h)^ {z) 



t^oo (tiQ — a)vQ 



(159) 
(160) 
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Limits of 



using 1 



lim e 

t-5-oo vq — a 



T^^ and S'(*): Setting Ak = dn (21)^/^ p, with Ak = 1, one finds, 
caling (117) and (118) for z, fc, m, with k—2m—l = (k— cj)/>(2t)^/^, 

^e-M.(t)T(i)(z-i) 



= lim ^_e-A^.(*) y % Afc 



t-5>oo — O 



fc>2m+l 



= lim 



t-5-oo — O 



fc>2m+l 



-k+2m 



i-vo) 



-k+2m 



-k+2m 

-fc+2m+l 



[(2t)i/3(i, 



fe>2m+l 

(161) 

using (153), the estimate for Qtin) and dominated convergence. In the same 
way 

/>oo 

hm {vo - a)e*^*Wr(2)(^) = -e^^'^'i / Q(K)e-'^2V3c^^ ^ -e2-fQ(-C). 

(162) 

Then for the 5^*^'s, one checks, using (157) and (159), the following limit: 

fc>2m+l\ ^0 a J 



lim 



fc>2mH 

/ Q(K)(iK / e-^'^''^^Ai{K + /3)(i/3 = y(C) 



and 



(163) 



/•OO /"OO 

lim 5(2)(z) = - / Q{K)dK / e2'^'f^Ai(/v + P)df3 = (164) 

To justify the limit, one uses the exponential estimate for Qt{K); also one 
uses an exponential estimate obtained from (159), (160) and ( |130 ) in order 
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to apply dominated convergence in the above integral. This ends the proof 
of Lemma 16. 6[ ■ 



'ext 



7 Scaling limits of the kernel 
7.1 Scaling limit as a perturbation of the Airy kernel 



In this section we prove the first part of the main Theorem 1.2, namely 



formula (19) 



Theorem 7.1 The following scaling limit holds under (137): 



lim {—VoY 



-ir-yK'^W^r,x;2s,y)pt'/' 



'-Tl>T2 



+ 2 



1/3 



e '1(ti-T2) 

Y^47r(Ti - T2) 



(165) 



where the first line on the right hand is the extended Airy kernel, as in {15), 



except for some conjugation; see (20). 



Proof: At first note that the y-scaling implies dy = pt^/^d^2- The kernel 
(— l)^~^K^^^(2r, x; 2s, y), as in (sj), is a sum of three parts. Multiplying each 
one with {—VQ)y~^~^^~^ pt^/'^ , the limit of the first part is taken care of by 



Lemma [6. 5| and the limit of the second part by Lemma 6.4 So, it remains to 
show the following uniform limit (uniform on Ci; ^2) ''"i) T2 in bounded sets). 



taking into account identity (152) and the symmetry of the Airy kernel: 



lim ptV3,(M--(i)-M-;i(i)) _ ^(2^^(0))^ia 



''—y,s J "—x,r 



P{2m+l,...) 



2I/3 



1 - X^KAiX^)-'Ag^^\ {X)All_^{X)dX, 



(166) 

remembering the notation ([s]) for the restriction of the kernel K{0) to [2m + 
1, 00). Using the identity 

(M-/^(i) -M-J^C,)) + {NUt)-M,{t)) = M^--^^^,,3.(i) - M--^^^,,3,(i), 
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together with Lemma 6.3 one checks the following limit: 



^Um ptV3,(M--(5)- _ i^«+i(0))-^a. 
2i/3 



?/,S! "—x,r 



£2(2m+l,...) 



17^ E 

f>2m+l 



(1 -e 



-M~^^ (-) 



2^/3 ^dA ((1 - X^ifAiX^)-iyi^7^^) (A)yig_^(A), 



where in the last equality, one uses the fact that 



Ai 



o(2t)i/3 



(167) 

dX, with Ae = l, 

and the fact, already pointed out just after (118), that i > 2m + 1 implies 
X>a = 22/3^. 

To justify the above, one needs uniformity of convergence. From the 



estimates (139), and using arguments, as in the derivation of the exponential 



bound for |Qt(K;)| in the proof of Lemma 6.6, we have that 



{product of the expressions in the first bracket in (167)} < cqc '^^ 



-ci A 



for S,i,Ti in a bounded set, and for A > a, with co,ci positive time- 



independent constants, while from (139) one has a similar estimate for the 



expression in the second bracket in (167). Thus by dominated convergence 



we have the convergence to the integral in (167), with uniform convergence 



for the parameters Ci,rj in bounded sets, concluding the proof of Theorem 



7.1, except for the identification with the Airy process kernel, which follows 
from (14). ■ 



7.2 Scaling limit as a sum of double integrals 



Here we indicate the proof of the second half of Theorem [L2l namely formula 
( [22I for r = s = n/2. To do so, we take the scaling limit of Kn^m{x, y), given 
by formula (102) in Proposition 5.2 for r = s = n/2;or what is the same, 
formula (113) containing L given by (112). It proceeds along similar lines 
as in |4J. We then merely sketch how to pass to the limit for the extended 
kernel ]K^^^(2r, x; 2s, y). 



68 



From (122), one has, using the scahng (117), ( 118| ) and (137) for for n, z, 

f$(z) = tF{z) + f/9(2t)i/3 log(-z) 

= iMs(t) + iC^-aC + 0(t-i/3)_ 



(168) 



With the scahng (118) and ( 137| ) for w,z,x,y, with Tj = 0, one finds the 
following limits: 



-w 



^y-l 



-z)y 1 



lim {-vof y I ' , + , , ^ , , ^ „ 



t^oo^ \ {-wY {-zY+^ 1 - - 



1 / e-^^ 



-I'D Ve 



+ 



and 



dzdw , 
lim ay 

t->oo z — 



-Wo- dr]. 

C-U! 



The £(x,y)-part of the kernel (113), 



-y) = L{x, y) + C(0; x - y), 



(169) 
(170) 

(171) 



consists of the four double integrals appearing in (112), with the Ei given 



in (100). Each of them will be examined separately, using the same saddle 



point arguments as given around Figure 18 in the proof of Lemma 6.2 



1st double integral: From (168) and (164), it follows that 



lim e5(*(^)-*("')) 



e3 



t— ^-oo 



1 s 

e3 



lim S'^^^z-^) = J'(C) and lim S^'^\w) = y{-uj) 



t— >oo 



and thus 



lim e^(*(^)-*("'))(l-5W(z-i))(l-5(2)(u;)) = ^i_(i_y(^))(i_y(_a;)). 



e3 



(172) 



2nd double integral: Since from (168), (163) and (164), 



lim e-*^-(*)e^(*(^)+*('")) 



n, JC'-C 



e3^ 
e 3 



lim {w - a)e^^*(*)T(2)(y;) = -e'^^^Qi- 
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one has 

lim (w - a)et(*(")+*("))(l-S(i)(z-i))r(2)(u;) 



e 3 



(i-nc))Q(-^). 

(173) 



3rd double integral: From (168), (162) and (164)), one obtains: 

lim e^*We-^(*(^)+*("')^ 



e3 



1 



lim 

t->oo z — a 



-2(tC 



Q(C). 



Therefore 



lim J_e-i(*(-)+*(-))r(i)(z-i)(l-5(2)(«;)) 
t^co z — a 



4th double integral: Since 



e 3 



1 3 



Q(c)(i-n-^))- 

(174) 



lim e 2 
t— >-oo 



($(tu)-<J>{z)) 



e3 



e3 



one has from (162) 



lim ^^e-t(*('")-*(^))T(i)(t.-i)r(2)(z) 
t-^oo a — w 



eh 

1 

e3^ 



-Q(-C)Q(u;). (175) 



Assembling the four double integrals, one finds the following limit, 
upon using the contours in [4J, in particular Figures 4 and 5: 

lim {-voy"" pt^^^ ^{x,y) 



/ du! — 3 

(5+iE i-<5+iE eT"~°"'^ 



(i-J'(C))(i-J'(-^)) 



dC / du! 

2(5+iK Js+iR e~°3" 



e 3 



+ 



f'X (1 - J'(C))Q(-c.) 



(271 



/ dw 



e 3 



+ 



(5+iR J-2(5+iK gT" 



f \ (1 - J'(-a;))Q(C) 



- dC duj 



e 3 



+ 



Q(-C)Q(^) 



(176) 
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This leads to the double integral representation (22) in Theorem 1.2, for 
the non-extended kernel r = s = n/2, upon interchanging uj and in the 
last integral. 

In order to take the scaling limit of the C{0;x — y)-part of the kernel 



(171), one needs to express C{0;x) in a more convenient form, as is given 
by formula (184) in Proposition 9.1 of the Appendix. Using this form, one 
now proves the following: 



Proposition 7.2 Under the scaling (131) for x,y with ti = T2 = and 
^1 I— 7- ^ and r]i i— )■ ij, one finds 

lim i-vof-ypt^/'^ C(0; x - y) 
= 2-1/3 

j dK Q{k) (^Ai(K + - r/)2-i/3) + Ai(K - (e - r])2-^/^)j 
+ dK d\ Q{k)Q{X) 

'Kj,,{k, a - - 7?)2-i/3) + K^^^i^ + (e - 7?)2-i/3^ A)^ 

Proof: Setting 1 = Ay = pt^^^dtj and 1 = Ak = pt^^^dn, one checks in the 
first expression (|184l), using (|129|) and (|153l), that: 



fc>2m+l 

^hm2-V3ePl°g(-o)(?-r?)<^/^p(2t)i/3rf^ ^ q'^^^ g'^^l^_^p{2t)'/^dK 

fc>2m+l 

>o™ I 1 V U / 



A,'>2m+1 



(177) 



= 2-^/^dr] / Q(K;)Ai(K + - 77)2-^/3)^ 
using dominated convergence in the integral. Similarly 

/"OO 

lim(-^;o)^-^Ay i^.+y = '^'^'''dn dK Q{K)Ai{K-{C-v)2-' 

fc>2m+l 



/3^ 
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The limit of the last two expressions in (184) are the following: 



k/>2m+l 

(1)^(2) j^(l) 

I —(70 J - t^ll 7 7 W. 

t—^OO 



k>2m-\-ll>2m+l 



poo poo 

/ dK i dX Q{k)Q{X)Km{k, A - - r])2'~^/^)dK dX 



and 

i(i)n(2)j^(i) 



k/>27n+l 

f'OO fOO 

= 2^1/3^7?/ dn dXQ{K)Q{X)K^,{K + {^-T])2^^/^,X), 

J U J (7 



thus establishing Proposition 7.2 



8 The tacnode kernel and non-colliding Brownian 
motions 



In this section, we prove expression (21 ) for the kernel K as given in (23), 



and we also show that K*^'^ = To do so, introduce the function 

= AiM(e-a + 2i/3^) 
and the operator T on L'^(a, oo), with kernel 

r(K,/3) = Ai(K + /3-a). 
Note that Kai, as an operator on L'^{a, oo), is given by 

/•oo 

/ Ai(K + A-CT)Ai(/3 + A-a)dA = i^Ai(K,/3), 

J a 

and that 

Km = T\ (178) 
It therefore follows that the resolvent of the Airy kernel can be expressed as 

oo 

1 + i? = (1 - KmY^ = (1 - r2)-i = T^^ (179) 

r=0 
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With this notation, expression (14) for A'^{k) can be written 

AI^^{k) = SI{k)-TSL^{k). 

Setting 

/■oo 

{f{K),g{K)) = / f{n)g{K)dK, 

J a 

one checks that 



(180) 



(181) 



and, using (179), 



1^ ((1 - X,KmX,)-UI[_^) {X)A^;^^{\)d\ 

oo oo 



oo oo 



(182) 







To conclude, by adding the two contributions (181) and (182), the kernel. 



as in ( 19 ) takes on the foUowing form: 

lfC*"^(Ti,e2;r2,6) 
= -'^Ti>T2P{n - -r2;6,6) 



^— T2 \ 

-6/ 





oo oo 


lri>T2P(n - T2;Cl,6) 

'(l-E:Ai)^i5[^)(A)S-;^(A) 



+ 2^3 / dA 



(l-i^TAO^W^^ J(A)5^--(A) 



(183) 



which yields formula (21), upon using formula (180) for ^^^_^(A) 
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Proof of Theorem 1.3: The kernel as in (23) can be written, using the 



notations (179) and (180) as follows 



[*br' = -lri>r2P(ri-T2;ei,6) 



+ 



2i/3 J 



1 + R){SZ^^-TS^P , TSll 
= -lTi>r2P(n - T2;6,6) 

oo oo 





oo 



which yields the formula in (183), ending the proof of Theorem 1.3 



9 Appendix: An expression for the integral C(0; x) 
in the kernel 



As was pointed out in Proposition 7.2, in order to take scaling limits, one 
needs to express C(0;x), defined in ( [lOTt , in a more convenient form. This 
is done below. In this section, let C{x) := C(0;x) and Cj(x) := Cj(0;x). 

Proposition 9.1 .■ The following holds: 

c(o;.)= E (ct^'sfi + er^g 

fc>2m+l 
k>2m+l 



(184) 



Proof: In this proof one uses over and over the expressions ( 72 ) for R^^^ {z ^ 



and R^'^\w), expressions for S'(^) and r(^), the series expansions m8h 
for h^^^ and the integral representation i52i for (7^*^ = gf\n). Using 



ro 



dz 
2TTi{-z) 
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grals for gf' , 



one checks by substituting the expressions above, and recognizing the inte- 



fc>2m+l fc>2m+l 



(185) 



and 

/■ dz 



k,£>2m+l kl>2m+l 
a,/3>0 ~ 



(186) 



k,e>2m + l 
a>0 



So, Ci(x) has the following form: 

Ci{x) = S^^oCUx) (187) 

with (adding (^^^=0 in the expression below is harmless, since CJ'(x) gets 
multiplied with 6x=/=o), 



k>2m.+l fc>2m+l 

^(1)^(2) (2) (1) . n(i)n(2)A 

fc,£>2m + l fc,^>2m+l 
a,/3>0 ~ 



-X 

k,t>2m + \ 
a>0 



(188) 



Also, by straightforward substitution, one finds: 

C2{x)= Yl Ql'^Q?^^/-k- (189) 

k,e>2m+l 

In the Lemma below, it will be shown that C^(0) = and thus Ci(x) = 
Cl{x). Therefore, adding (188) and (189), one checks, using (54) and the 
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identity ( 192 ) of the Lemma below: 



Cix) = Ciix) + 2C2(x) = Clix) + 2C2(x) 

^ A:>2m+1 ^ ^ k>2m+l 

fc,£>2m+l 

\^ ^(1)^(2) (2) (1) . v;^ n(^)n(2)A 

fc,<>2m+l fc,£>2m+l 

a,;3>0 

fc>2m+l 



>2m+l 



fc>2m+l 



(190) 



This estabUshes Proposition 9.1 



Lemma 9.2 The following holds 

C7*(0) = (191) 

and 

*:,«>2m + l fc/>2m+l 

fc>2m+l fc>2m+l 

(192) 

Proof: (i) First we prove that 

E off!" = E 

k>M k>M 

Since, by definition, upon setting := (1 - XmK'-^\0)Xm)~^ , (set 2m + 
1 = M) 

= ((1 -XMi^(2)(0)XM)-^W)) (fe) =: P(')(W^)(A:) 
= ({l-XMK<^'H0)XMr'XMg^^A (k) =: P(i)(Xm5^'^)(A:), 
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one has, using P^^^~^ = p(^), 

k>M 



k>M 



{ii) Therefore, by ( |188[ ), 



(i)„(2) , >^ ^(l)r1(2)^^(l)^n^ >^ ^(1)^(2) 

k 

k>M kl>M k>M 



+ {P^^H^M9^'^),^MK^'H0)XMP^'HXMg^^^)) (194) 



Then CJ'(O) = impUes Ci{x) = CJ'(x), estabUshing identity (191). 



(iii) Using ' , for k > M, 

Qi"= 9l" + E4?(0)<3j" 

£>M 

Qf'= 9f+E<(0)<3f 

£>M 

and writing 1 = {5x>o + ^Sx=o) + {Sx<o + ^Sx=o), 

A:,£>M 

= ^ Qfc^^Q£^'''5x,^-fc (('^x>0 + 5'^2:=o) + {Sx<0 + 55x=o)) 
fc,f>Af 



(195) 



'for P = (1 - K)-'^, one has P = 1 + ifP 
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(2) 

I 



\k>M k,i>M 
,k>M kl>M 



(196) 



Also, by (54), 



k,i>2m + l 



(197) 



kl>M 



Hence, subtracting these two formulas, one finds 

\^ ^(1)^(2). v;^ ^(1)^(2) (2) (1) . 



kl>M 



k,l>M 
a,l3>0 



(2)„(1) 



(198) 



A;>M 



thus ending the proof of identity (192) and Lemma 9.2 
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